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Abstract. In this paper we formulate and analyze the Multi-Configuration 
Time-Dependent Hartree-Fock (MCTDHF) equations for molecular systems 
with pairwise interaction. This is an approximation of the TV-particle time- 
dependent Schrodingcr equation which involves (time-dependent) linear com- 
bination of (time-dependent) Slater determinants. The mono-electronic wave- 
functions satisfy nonlinear Schrodinger-type equations coupled to a linear sys- 
tem of ordinary differential equations equations for the expansion coefficients. 
The invertibility of the one-body density matrix (full-rank hypothesis) plays 
a crucial role in the analysis. Under the full-rank assumption a fiber bundle 
structure shows up and produces unitary equivalence between different useful 
representations of the approximation. We establish existence and uniqueness 
of maximal solutions to the Cauchy problem in the energy space as long as 
the density matrix is not singular for a large class of interactions (including 
Coulomb potential). A sufficient condition in terms of the energy of the ini- 
tial data ensuring the global-in-time invertibility is provided (first result in 
this direction). Regularizing the density matrix breaks down energy conserva- 
tion. However a global well-posedness for this system in L 2 is obtained with 
Strichartz estimates. Eventually solutions to this regularized system are shown 
to converge to the original one on the time interval when the density matrix 
is invcrtible. 



Multi-configuration methods, Hartree-Fock equations, Dirac-Frenkel variational 
principle, Strichartz estimates 



The purpose of the present paper is to lay out the mathematical analysis of 
the multi- configuration time-dependent Hartree-Fock (MCTDHF) approximation 
which is used in quantum chemistry for the dynamics of few electron problems, or 
the interaction of an atom with a strong short laser-pulse [7J [37J [3H] and [2T] ■ The 
MCTDHF models are natural generalizations of the time-dependent Hartree-Fock 
(TDHF) approximation, yielding a hierarchy of models that, in principle, should 
converge to the exact model. 

The physical motivation is a molecular quantum system composed of a finite 
number M of fixed nuclei of masses mi , . . . , mjj > with charge z\ , . . . , zm > 
and a finite number N of electrons. Using atomic units, the iV-body Hamiltonian 
of the electronic system submitted to the external potential due to the nuclei is 
then the self-adjoint operator 



acting on the Hubert space L 2 (£l N ;C) with pairwise interaction between the elec- 
trons of the form 



1. Introduction 



(1.1) 






l<i<j<N 
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with v real- valued and v > 0. Here and below f2 is either the whole space R 3 
or a bounded domain in R 3 with boundary conditions. The N electrons state is 
defined by a wave-function ^ = ^(xi, . . . ,xn) in L 2 (fl N ) that is normalized by 
H^Hz^fiiV) = 1. To account for the Pauli exclusion principle which features the 
fermionic nature of the electrons, the antisymmetry condition 

^(x lt x N ) = e(o-)*(x (r (i) , . . . , x^ N )), 

for every permutation a of {1, ... , N} is imposed to the wave- function The 
space of antisymmetric wave- functions will be denoted by Ai=i L 2 (£l). In and 
throughout the paper, the subscript n of — A Xi means derivation with respect to 
the i th variable of the function \t. Next, 

M 

U(x) := - V 



^— ' \x - R m \ 

m— 1 1 1 

is the Coulomb potential created by M nuclei of respective charge Zi, • ■ ■ , zm > 
located at points - ,Rm €E M. 3 and v(x) ~ is the Coulomb repulsive 

potential between the electrons. Actually our whole analysis carries through to 
more general hamiltonians (possibly time-dependent) as explained in Section [7] 
below. 

For nearly all applications, even with two interacting electrons the numerical 
treatment of the time-dependent Schrodingcr equation (TDSE) 

(1.2) i?l=H N y, tt(0) = *°, 

at 

is out of the reach of even the most powerful computers, and approximations are 
needed. 

Simplest elements of AiLi L 2 (Sl) are the so-called Slater determinants 

(1.3) V(xi,...,x N ) = -=det((f»i(x j )) lt 



/JV! V JJ Jl<i,j<N 

constructed with any orthonormal family fa in L 2 (il) . The factor ensures the 
normalization condition on the wave-function. Such a Slater determinant will be 
denoted by fa A . . . A 4>n ■ The family of all Slater determinants built from a complete 
orthonormal set of L 2 (VL) is a complete orthonormal set of /\ i=1 L 2 (Q). Algorithms 
based on the restriction to a single Slater determinant arc called Hartree-Fock 
approximation (HF). On the other hand the basic idea of the multi-configuration 
methods is to use a finite linear combinations of such determinants constructed 
from a family of K(> N) orthonormal mono-electronic wave- functions. 

One observes (this computation is done in Subsection [53]) that in the absence of 
pairwise interacting potentials any Slater determinant constructed with orthonor- 
mal solutions fa(x,t) to the single-particle time-dependent Schrodinger equation 
gives an exact solution of the A-particle non interacting time-dependent Schrodinger 
equation. Such fa(x,t) are called orbitals in the Chemistry literature. The same is 
true for any linear combination of Slater determinants with constant coefficients. 
Of course, the situation turns out to be completely different when pairwise inter- 
actions are added : a solution to TDSE starting with an initial data composed 
of one or a finite number of Slater determinants will not remain so for any time 
t 7^ 0. Such behavior (called "explosion of rank") is part of the common belief, but 
is not shown rigorously as a property of the equations, to the best of our knowl- 
edge. In the MCTDHF approach one introduces time-dependent coefficients and 
time-dependent orbitals to take into account pairwise interactions and to preserve 
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the finite linear combination structure of Slater determinants in time. Using time- 
independent orbitals as it corresponds to a Galerkin-type approximation would save 
the effort for the nonlinear equations, but requires a much larger number of relevant 
orbitals and hence the numerical cost is much higher. The motion of the electrons 
in the MCTDHF framework is then governed by a coupled system of K nonlinear 
partial differential equations for the orbitals and (^) ordinary differential equations 
for the expansion coefficients (see for instance System (|3.27|) ). 

Although MCTDHF is known for decades, the mathematical analysis has been 
tackled only recently. For a mathematical theory of the use of the time-independent 
multi- configuration Hartree-Fock (MCHF) ansatz in the computation of so-called 
ground- and bound states we refer to |25[ I17[ I26| . A preliminary contribution 
was given by Lubich [28] and Koch and Lubich [24) for the time-dependent multi- 
configuration Hartrec (MCTDH) equations for bosons, for the simplified case of a 
regular and bounded interaction potential v between the electrons and a Hamilton- 
ian without exterior potential U. The MCTDH equations are similar to MCTDHF 
from the functional analysis point of view, although more complicated from the 
algebraic point of view, since more density-matrices have to be considered in the 
absence of a priori antisymmetry requirements on the A^-particle wave- function (see 
also [23] for an extension to MCTDHF equations). Using a full-rank (i.e. invertibil- 
ity) assumption on the one-body density matrices, the authors proved short-time 
existence and uniqueness of solutions in the functions space ff 2 (R 3 ) for the orbitals 
with the help of Lie commutators techniques. Numerical algorithms are also pro- 
posed and analyzed by the groups around Scrinzi (e.g. |37| ) and Lubich, the proof 
of their convergence generally requires the H 2 -type regularity assumptions (see e.g. 

We present here well-posedness results for the MCTDHF Cauchy problem in H 1 , 
H 2 and L 2 , under the full-rank assumption on the first-order density matrix and for 
the physically most relevant and mathematically most demanding case of Coulomb 
interaction. We also give sufficient conditions for global-in-time full-rank in terms 
of the energy of the initial data. Eventually solutions to a perturbed system with 
regularized density matrix are shown to converge to the original one on the time 
interval when the density matrix is invertiblc. 

This paper is organized as follows. In Section 2 we give a complete analysis 
of the ansatz associated to the multi-configuration Hartrcc-Fock approximation. 
Essentially, this ansatz corresponds to a linear combination of Slater determinants 
built from a vector of complex coefficients C and a set of orthonormal, square 
intcgrablc functions represented by a vector $ = (<f>\, cf>2, . . . , 4>k), for K > N . The 
first-order density matrix is introduced and represented by a complex-valued matrix 
F which corresponds to the representation of the kernel of the first-order density 
matrix in the orthonormal basis {4>i, 4>2, ■ ■ ■ , 4>k}- By abuse of language this matrix 
T depending only on the expansion coefficients C is also called density matrix. 
Its invertibility is a crucial hypothesis which will be referred to as the full-rank 
hypothesis. Under this hypothesis, the corresponding set of pairs (C, $) is endowed 
with a structure of a fiber bundle. In Section 3, two set of equivalent systems are 
presented. The first one, iSo, called variational system is inspired by a variational 
principle. The second one, 5h, will be referred to as working equations. In Section 
4, the system Sq is used to prove the propagation of the normalization constraints, 
the conservation of the total energy and an a posteriori error estimate for smooth 
solutions (if they exist). The system Su is used to prove local existence, uniqueness 
and stability with initial data in H m for to > 1. In particular the space H 1 is used 
to balance the singularity of the potentials (of Coulomb type) and we prove the 
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local well-posedness using the Duhamel formula. Next, the conservation of the 
total energy allows to extend the local-in-time solution until the associated density 
matrix IT becomes singular. Therefore, Section 5 is devoted to a criteria based 
on the conservation of the energy that guarantees the global-in-time invertibility 
of the density matrix IT. To handle the possible degeneracy of this matrix, a 
regularized problem is considered in Section 6. For this problem the conservation 
of the energy does not hold anymore. Hence, we propose an alternate proof, also 
valid for singular potentials, but that is only based on mass conservation. Such 
proof relics on Strichartz estimates. Eventually, one expects that the solution of 
the regularized problem converges towards the solution of the original one as long 
as the unperturbed density matrix is invertible. The proof is a H 1 version of the 
classical "shadowing lemma" for ordinary differential equations. Finally, in the last 
section we list some extensions to time-dependent Hamiltonian including a laser 
field and/or a time-dependent external potential. The case of discrete systems is 
also discussed there. 

Some of the results presented here have been announced in [31] and [3] and the 
details of the L 2 theory are worked out in [ST] , 

Notation. (-, •) and (-|-) respectively denote the usual scalar products in L 2 (il) 
and in L 2 (fl N ), (•, •) the scalar product in L 2 (Q) K and a ■ b the complex scalar 
product of two complex vectors a and b. The bar denotes complex conjugation. We 
set L^(£l N ) := Afcli L 2 (Q) where the symbol A denotes the skew-symmetric ten- 
sorial product. Throughout the paper bold face letters correspond to one-particle 
operators on L 2 (fl), calligraphic bold face letters to operators on L 2 (£l N ), whereas 
"black board" bold face letters are reserved to matrices. C(E; F) denotes the set 
of continuous linear applications from E to F (as usual C(E) = C(E; E)). 
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2. Fiber Bundle Structure of the Multi-Configuration 
Hartree-Fock Ansatz 

2.1. The MCHF ansatz. For positive integers N < K, let Sjv.k denote the set of 
increasing mappings a : {1, . . . , N} — ► {1, . . . , K} 



■>N,K = {<T = {a(l)<...<a(N)}c{l,...,K}}, 



N 



For simplicity the same notation is used for the mapping a and its range {<r(l) < 
. . . < a{N)}. Next we define 



with 



(2.1) O iW = {* = Wi,..,fe)ei 2 (^ : J kfadx = 8 iJ }, 

with Sij being the Kronecker delta and with S r being the unit sphere in C r 
endowed with the complex euclidean distance 

(2-2) fiT- 1 = {C = {c a )^ N>K £C r : ||C|| 2 = J2 M 2 = l} 

a 

with the shorthand for X^es^ k' ^ an ^ a e ^n,k and $ in C^ 2 (n) K i we 
associate the Slater determinant 



<Mxi, . . . ,Xjv) = <Ar(l) A ... A CT (JV) 



AM 



4>cx(l){x\) ... 4>ct(1)(xn) 
4><r(N)(xi) ■■■ 4>a-(N)(x N ) 



The mapping 

(2.3) (C, $) .— > * = 7r W)J r(C, $) = ^ c ff $ c 



is multilinear, continuous and even infinitely diffcrcntiablc from J~n,k equipped 
with the natural topology of C x L 2 (fl) K into L\ (0^). Its range is denoted by 



A' = k(Fn,k) ={* = ]Tc a $ CT : (C,$) e Fn,k}- 



When there is no ambiguity, we simply denote 7r = -kn,k- The set Bn,n is the 
set of single determinants or Hartree-Fock states. Of course Bn,k C Bn,k' when 
if' > if and actually 

lim B w ,* = ftf eL^ft") : ||*|| = l}, 

in the sense of an increasing sequence of sets, since Slater determinants form an 
Hilbert basis of L 2 (fi^) (see [27])- In particular, for cr, r € Sat.k, we have 

(2.4) ($ CT I * r ) = cT CT , T . 
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Observe that without orthonormality condition the formula (|2.4p becomes 

(2.5) {fa A...A<p N | a A...Aftr) = det((^;e i )) 1 < iJ < JV 

for $, 3 G L 2 (r2) W which will be used below (see [2"7]). 

The set of multi-configuration ansatz Bn,k is characterized in Proposition 12. 21 in 
Subsection 12 . 21 in terms of the so-called first-order density matrix, and its geometric 
structure is analyzed in Subsection 12.31 

2.2. Density Operators. For n = 1, . . . , N and for * G L 2 A (Q N ) with ||*|| = 1. 
a trace-class self- adjoint operator [vP ® , called n th order density operator, is 
defined on L\ (Q n ) through its kernel ® \&] . 

(2.6) [* ® *] JX n , Y n ) = ( N ) [ (JT nj z^) W(y„, Z,f) dZ», 
for 1 < n < N — 1 and 

[* ® *] :JV (Xjv, Fjv) = *(Xjv) *(yjv), 

with the notation 

-A n (^1 j ■ ■ ■ 3 ^n)) (l/l 3 * • • j yn)} ; 

Z„ = (Zn+lt ■ ■ ■ 7 z w), dZ n = dz n+ i . . . dzpj, 

and similarly for other capital letters. Our normalization follows Lowdin's |27j . A 
simple calculation shows that, for 1 < n < N — 1, 

(2.7) [*®*] :n (x nj y„) = -£±-i y [*<g>*] :n+1 (x n ,*,y n) 

In particular, given l<n<p<N— 1, one can deduce the expression of [\f ® \f] . 
from the one of [\& (8) . These operators satisfy: 

Proposition 2.1 ( fTJ [T31 HH [22] ) • For every integer 1 < n < N, the n-th order 
density matrix is a trace-class self-adjoint operator on L\(£l n ) such that 

(2.8) 0<[*®¥] in <l, 
in the sense of operators, and 

Tr iW [*®*] !B = 

Actually, multi-configuration ansatz correspond to first-order density matrices 
with finite rank, and we have the following 

Proposition 2.2. [Lowdin's expansion theorem [57]; see also [TTJdS]/ Let K > N, 

then 

B N , K = Tr(F NtK ) = {V € L 2 A (n N ) : ||#|| = 1 and rank[* ® < K}. 

More precisely, if * = 7r(C, $) with (C, $) € ^~n,Kj then rankf* ® *] < K and 
Ran[* ® *] C Span{0i, • • • , (p K {- If * G Sjv.k and if rank[* ® \&] = AT' with 
N < K' < K and with {0i, . . . , <^r-'} being an orthonormal basis of Ran[^ ® 
then \1/ can be expanded as a linear combination of Slater determinants built from 
{(pi, - • • , <Pk'}- The first-order (or one-particle) density matrix [\& <g> \I>] is often 
denoted by 7^ in the literature and in the course of this paper. According to 
Proposition ^ . II above it is a non-negative sclf-adjoint trace-class operator on L 2 (Vl), 
with trace N and with operator norm less or equal to 1. Therefore its sequence 
of eigenvalues {7i}i>i satisfies < 7, < 1, for all i > 1, and X)j>iTt = N. In 
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particular, at least TV of the 7i's are not zero, and therefore rank 7$ > TV, for any 
* G L*(Sl N ). 

Similarly, if * = tt(C, $) G Bn,k, the range of the operator [ir(C, $) <g> ir(C, $)] : „ 
is n Span{$} and its kernel is ^ ® n Span{$}^ with Span{$} := Span{^i, . . . , (/>#}. 

Therefore the operator is represented by an Hcrmitian matrix in n Span{<£>} whose 
entries turn out to depend only on the coefficients C and the dependence is qua- 
dratic. For the first- and second- order density operators we have the explicit 
expressions 

Proposition 2.3 ([17], Appendix 1). Let ^ = 7r(C, $) in Bn,k, then the operator 
kernel of the second-order density matrix kernel is given by 

K 

(2.9) [V®*] :2 (x,y,x',y') = ^ 7yw &0=) <t>M ? fc (a/) W) 

i,j,k,l=l 

with 



(2.10) ^ = 1(1-^X1-4,0 Yl (-i)^(-i)Lc CT c T , 

where for i ^ j. 



<T,T I i,j£(T, fc,iGT 
<?\{i,j}=T\{k,l] 



Similarly, the kernel of the first-order density matrix is given by the formula 

K 

[*® *]..i(x,y) = ^ in M*)4j(v) 



with 



(2.11) ™ = ^f>w* = Yl (-i) CT " 1(i)+r " Iw c^ 

fc — 1 (T,T I i£(7,j€iT 

<r\{i}=r\{j} 

and 

(2.12) 7 «= ^ |c,| 2 . 

cr I i£er 

The first-order density matrix allows to characterize the set Bm,k (see Propo- 
sition 12.21 above) whereas the second-order density matrix is needed to express 
expectation values of the energy Hamiltonian as soon as two-body interactions are 
involved. 

Since the coefficients 7^ only depend on C, we denote by F(C) the K x K 
Hermitian matrix with entries 7^, 1 < i,j < K (the adjoint of the matrix rep- 
resentation of the first-order density operator in Span{$}). The matrix F(C) is 
positive, hermitian and of trace TV with same eigenvalues as 7^ and same rank, and 
there exists a unitary K x K matrix U such that U T(C) U* = diag(7i, . . . ,7k) 
with < -fk < 1 and J2k=i Ik = N. Hence, 7* can also be expanded as follows 

K 

(2-13) l9 {x,y) = Y,!i4>'i{x)Wi{y), 

where = C/ • $ with obvious notation and with {(/)[, • • • , being an eigenbasis 
of 7$. Note that it is easily recovered from (|2.12[) that < 7, < 1 for C G S' r_1 . 
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Remark 2.4. When K = TV (Hartree-Fock case), 7$ being of trace TV must be the 
projector on Span{$}; that is 

N 
i=l 

with P$ denoting the projector on Span{$}. In this case, (|2.11[) and (|2.10[) simply 
reduce to 7^ = S ltj and -f ljk i = %(5i t k8jj - 5 it iSj^); that is T(C) = In- 

The representation of a wave-function \& £ Bn,k in terms of expansion co- 
efficients C and orbitals $ is obviously not unique as it is already seen on the 
Hartree-Fock ansatz. Indeed, if ^ HF = <f) 1 A • • • A 4>n — ipi A • • • Ai/jn, there exists a 
unique TV x N unitary transform U such that (4>i, • • • , 4>n) = (ipi, • • • , iPn) • U . The 
pre-image of ty HF by 7r in Tn,n is the orbit of (<pi, ■ ■ ■ ,(/>n) under the action of 
On, with Oe being the set of £ x t unitary matrices. In the general case under the 
full-rank assumption the set Bn,k has a similar orbit-like structure as explained 
now. 

2.3. Full-rank and fibration. We introduce 

9Bn,k '■= £ Bn,k ■ rank 7^ = A'} 

and, by analogy, 

dT N ,i< = 7T N ^ K (dB N , K ) := {(C, $) G T n ,k ■ rankr(C) = K}. 

dJ^N,K is the open subset of J-n,k corresponding to invcrtiblc F(C)'s (full-rank 
assumption) . 

Clearly 8Bn,n = Bn.n and OTn.n — Fn,n\ that is, the full-rank assumption is 
automatically satisfied in the Hartree-Fock setting (see Remark l2.4p . 

On the opposite, it may happen that 8Bn,k = (in that case Bn,k = Bn,k-i)- 
Indeed, for K > N the admissible ranks of first-order density matrices must satisfy 
the relations Q71 H>] 

f = 1 TV = 1 

K < > 2, even TV = 2 . 

I > TV,^TV+1, TV > 3. 

From now on, we only deal with pairs (TV, K) with A admissible. We recall from 
[2"T] the following 

Proposition 2.5. Let (C, $) and (C, $') in 8T n ,k such that tt(C, $) = tt(C", $')• 
Then, there exists a unique unitary matrix U £ Ok ond a unique unitary matrix 
d(U) = U £ O r such that 

$' = [/•$, C = d(U) ■ C 

where, for every a £ Sjv,ifj 

T 

Moreover, 

(2.14) T(C) = UT(C) U*. 

Proof. Let (C,$) and (C",$') in dT N ,K such that tt(C, $) = tt(C*',$') = * £ 
8Bn,k- From Proposition 12. 2[ Span{$} = Span{$'} = Ran(7*) with $ and $' 
in Ojjuq^k. Therefore, there exists a unique unitary matrix U £ Ok such that 
$' = [/ • Eqn. (|2.14p follows by definition of F(C). Accordingly, there exists a 
unique unitary matrix U in O r that maps the family {$ T } rS s N K to {^' a } ae ^ N K . 
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More precisely, being given a G ^n,k, we have by a direct calculation (see also 



(2.15) 



where, for all er, r £ Ej/ t K, 

^t(1),t(1) 



(2.16) 



a(JV),r(l) 



= det(C7, 



"'0')>''"(i)Ji<. I -j<jv 



U a (l),T(N) • ■ ■ £/cr(jV),T(j\T) 

= det((0' CT(j . ) ;0 T(i) )) i ^^ jv . 

By construction the rxr matrix U with matrix elements U CTjT is unitary. By the 
orthonormality of the determinants, we have 

(2-17) c' a = (ir(C, $) | <) = J2 c t(®t I K) = Cr, 

r r 

whence the lemma with rf(?7) = U. □ 

Under the full-rank assumption and given (N, K) admissible, the set 8Bn,k 
is a principal fiber bundle. In differential geometry terminology, OBn.k is called 
the base, and, for any '5 € dBjy,K, the pre-image 7r — 1 (^) is i/ie fiber over "J. 
Proposition 12 . 51 defines a transitive group action on OJ^n^ according to 

(2.18) 



(C, $') = U • (C, $) <^=> C" = d(E/) • C and $' = [/•$, 
W := (d(U),U) eO r x K - 



Indeed on the one hand, it is clear from the expression for the matrix elements of 
d(U) that dfljt) = I-r- On the other hand from (|2.15[) and (|2.17p it is easily checked 
that d(UV) = d{U) d(V). Therefore couples of the form (d(U), U) form a subgroup 
of O r x Ok that we denote by O r K . The action of O r K is not free on Tn,k itself — 
this is illustrated in Remark 12.61 below on the examples of Slater determinants in 
Fn,k with K > N — , but it is free on 8Tn,k and transitive over any fiber 7r _1 (^') 
for every ^ G OBn.k- Therefore, the mapping tt defines a principal bundle with 
fiber given by the group O r K . We can define local (cross-)sections as continuous 
maps s : * i— ► (C, $) from 3Bn,k to 8Tn^k such that tt o s is the identity. In 
particular, dTN,K/O r K is homeomorphic to 8Bn,k- Since the map tt is C°°, one 
concludes from the inverse mapping theorem that the above isomorphism is also 
topological. In the Hartree-Fock case K = N where the full-rank assumption is 
automatically fulfilled, tt^j 1 n {Bn,n) is a so-called Stiefel manifold. 

Remark 2.6. The following example illustrates the necessity of the full-rank as- 
sumption. As 

K<K' B N>K CB N , K ', 

any Slater determinant fy HF ~ (j>i A • • • A 4>m can also be seen as an element of 
Bn,k for all K > N . If K > N , the pre-image of fy HF by tt in Tn,k does 
not have a similar orbit structure as shown on the following example. Let C = 
(1,0, • . • , 0) G S r where all coordinates but the first one are and let $' = 
(0i, . . . ,(j> N ,4> N+1 , . . .,0k) e L 2 in) K with (f>i G Span{0i, . . . , (/W} 1 - for every N + 
1 < i < K, then (C, $') G Tn,k and ^hf = tt(C, $'). There is no group-orbit 

structure on I Span{0i, . . . , <pN 



K-N 
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Having equipped 8Bn,k with a manifold structure we turn to the study of its 
the tangent space. 

Being multi-linear with respect to the variables C and $, the application 7r is 
clearly infinitely differentiable. Its gradients 

Vtt: C r x L 2 {Q) k — ► C(C r ;L 2 (n N )) x C(L 2 (n) K ;L 2 (Q N )) 

* = 7r(C,$) i — ► V = (V c V* *) 

are computed as follows for every (C, $) G Tm.k '■ 

i) for any <5C in C r , 

r ~ j r 

(2.19) VcnSC]=J2^k^— = J2 Sc ^^ 

k=l 0cak fe=l 

with Sjv.k: = {<ti, • • • ,oy}, 

ii) foranyC=(Ci,---,C*:)e£ 2 (n)* 

(2.20) v^[ C ]=E#[c fe ]= E C *£|^ 

with 

9* r 

(2.21) — [<,•] =V(iW / ¥(x 1 ,...,x N )(l> k (xi)dx i . 
0(pk i=l Jn 

Remark 2.7. For every a G ~^n,k and every 1 < k < K , we have 
(2.22) 

<9$<r f A ' ' ' A ^(i-l) A C A CT (j+l) A • • • A ff(A r) «/ cr" 1 (fc) = j, 

if kg a. 



■[(} 



Remark 2.8. From (|2.22p we recover the Euler Formula for homogeneous func- 
tions, that reads here 

(2.23) *=lf;^]:=Iv.*[*]. 

From the definition of the adjoint V$** G £(L 2 (f^); L 2 ^)^) of the operator 
V$<I' one has 

(2.24) VC G L 2 (Qf, VS g L 2 A (Q N ), (V***[S] ; 0z»<n)* = (s| V** [C]) , 



with 



— — [E](x) = N (f> k (y)( E(x,X2,...,x N )'$r(y,X2,... > x N )dx2---dx N )dy 
o<Pk Jn / 

for all 1 < fc < J\, for any function S in L 2 (il N ). 

It is also worth emphasizing the fact that changing (C, <j>) to (C, $') following the 
group action (|2.18p , involves a straightforward change of "variable" in the derivation 
of Sb; namely, with a straightforward chain rule, 

(2.25) V c # = U* • Vc * = Vc* • d{U), V$ * = V$/ * • C/ 
and 

(2.26) [V$/ *]* = ?/• [V$ *]*. 
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The following further properties of the functional derivatives of "J will help to link 
the full-rank assumption with the possibility for tt to be a local diffeomorphism in 
a neighborhood of = k{Cq, <&o) G 9Bn,k- 

Lemma 2.9. Let (C, $) G F N , K with f = n(C, $). Then, for all ( G Spanf*}- 1 , 
£ G L 2 {VL) and cr, r G £;v,K", we /lave 

and 

/or any 1 < fc, / < A'. 

Proof. The first claim follows immediately in virtue of (|2.22|) and (|2.5p . For the 
second claim we proceed as follows. Thanks to (|2.22p again 

a* a* \ ^ _ /9$ CT 9d> T 

2 (-ir 1 «(-ir 1(i) c CT c r (c ! o 

<7,T | fc£er, Z£r 
<A{k}=r\{l} 



We conclude with the help of IpTTj) , □ 

From l[2"25j> . ([TT^j) . jgHQ]) and (|2~2"Tj) . the tangent space ofdB N>K at * = tt(C, $) 
is given by 

(2.29) 1 CT CT fc=i 0fc 

<SC = (<5c CTl , • • • , 8c ar ) G C r , G Span^}- 1 , for every l<k<K 

Note that the tangent space (|2.29p only depends on the basis point ^ and not 
on the choice of coordinates (C, <&) in the corresponding fiber. In Physicists' ter- 
minology this is the space of allowed variations around (C, $) in J~n,k according 
to the constraints (|2.ip and (|2.2p on the expansion coefficients and the orbitals 
respectively. 

Let \&o = 7r(Co,$o) be in Bn,k with invertible T(Cq). Then the local mapping 
theorem at (Co,$o) allows to define a so-called section ir^ 1 : ^ i— > (C, <f>) as a C 1 
diffeomorphism in the neighborhood of ^q. According to (|2.29p . we have to check 
that (0, 0) is the only solution in C x Spanj^o} -1 to 

(2.30) dn (Co ,* o) (dC,S$) = Y,**Sc« + ^ EJz^]= - 

cr k— 1 

Indeed, on the one hand, if we scalar product the above equation with <E> r for any 
r G Stv.k we obtain o*C = in virtue of the orthonormality of Slater determinants 
and (|2.27p . On the other hand, for a given 1 < I < K and any £ G L 2 (f2), the 
scalar product of (|2.30p with yields 

E<H ^ I H ra) =|> - ((™)«'*> = 
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thanks to (|2.28[) . Since £ is arbitrary in L 2 and since T is invertible this is equiv- 
alent to 5$ = 0, hence the result. The full-rank property is mandatory for lifting 
continuous paths 1 1— > ^(t) on the basis 8Bn,k to continuous paths 1 1— * (C(t), $(£)) 
on 8Tn,k- 

2.4. Interpretation in terms of quantum physics. The wave-function \I> g 
L 2 (fi A ') with H^ll = 1 is interpreted through the square of its modulus ^(X/v)! 2 
(= [\& (g) \&] (Xn,Xn)) that represents the density of probability of presence of 
the JV electrons in Q N . More generally, for all 1 < n < N, the positive function 
X n ® (JT„,X„) is in L^Sl") with L 1 norm equal to (^), and it is 

interpreted as (^) times the density of probability for finding n electrons located 
at X n G 0". Any set {c(l), . . . , cx(JV)} for a S £jv,.k- is called a configuration in 
quantum chemistry literature and this is where the terminology multi- configuration 
comes from for wave-functions in Bm,k- When {4>k}i<k<K is an orthonormal basis 
of Ran[\I'(£>\E r ] each mono-electronic function <pk is called an orbital of 'P. When the 
orbitals are also eigenfunctions of [tv(C, (8) 7r(C, $)] : i according to (|2.13p they are 
referred to as natural orbitals in the literature whereas the associated eigenvalues 
{li}i<i<K are referred to as occupation numbers. Under the full-rank assumption, 
only occupied orbitals are taken into account. The functions with JV — 1 variables 
^{xi, . . . , xis[)<j) k {xi)dxi that appear in formula (|2.2ip are known as a single-hole 
function (see e.g. [SJ E] ) ■ Finally, the K x K matrix T(C) is called the charge- and 
bond matrix (see Lowdin |27jV 

"Correlation" is a key concept for many-particle systems. Whereas the "cor- 
relation energy" of a many particle wave-function associated to a many particle 
Hamiltonian is a relatively well-defined concept, the intrinsic correlation of a many 
particle wave- function as such is a rather vague concept, with several different def- 
initions in the literature (see among others [2"01 fTT)] and the references therein). 
In [19] Gottlieb and Mauser recently introduced a new measure for the correla- 
tion. This non-freeness is an entropy-type functional depending only on the density 
operator [\& ® ^]a, and defined as follows 



(£(*) = -Tt|[*®*]. 1 log([*®*]:i)|-Tr|[(l-[*®*] : i)log(l-[*®*]:i; 

Hence it depends on the eigenvalues of g) in the following explicit way 



= - E (n + (! - 7.) log(l - 7i)) 

j=i ^ ' 

It is a concave functional minimized for 74 = or 1. In the MCHF case this 
functional depends implicitly on K and N via the dependency on the 7^. This 
definition of correlation has the basic property that the correlation vanishes if and 
only if ^ is a single Slater determinant. The simple proof is based on the Lowdin 
expansion theorem (see Proposition 12.21 and Remark 12. 4p . 

The single Slater determinant case is usually taken as the definition of uncor- 
rected wave-functions. The Hartree-Fock ansatz is not able to catch "correlation 
effects" . When there is no binary interaction the Schrodinger equation propagates 
Slater determinant (see Subsection 13 . 5|) . However, the interaction of the particles 
would immediately create "correlations" in the time evolution even if the initial data 
is a single Slater determinant, - however, the TDHF method forces the dynamics 
to stay on a manifold where correlation is always zero. 

Improving the approximation systematically by adding determinants brings in 
correlation into the multi- configuration ansatz. Now correlation effects of the many 



MCTDHF EQUATIONS 



13 



particle wave- function can be included in the initial data and the effects of dynami- 
cal "correlation - decorrelation" can be caught in the time evolution. This is a very 
important conceptually advantage of MCTDHF for the modeling and simulation of 
correlated few electron systems. Such systems, for example in "photonics" where 
an atom interacting with an intense laser is measured on the femto- or atto-second 
scale, are increasingly studied and have given a boost to MCTDHF (see e.g. [7], [2]). 



3. Flow on the Fiber Bundle 

In this section, we consider a general self-adjoint operator Ti in L 2 (Q, N ). Most 
calculations here stay at the formal level with no consideration of functional anal- 
ysis. Solutions are meant in the classical sense and in the domain of the operator 
Ti . In Section H] below physical problems will be considered and details concerning 
proof of existence, uniqueness of solutions and blow-up alternatives in the appro- 
priate functional spaces will be given. 

From this point onward, T > is fixed. A key point of the time-dependent case 
is that the set of ansatz Bn.k is not invariant by the Schrodingcr dynamics. It is 
even expected (but so far not proved to our knowledge) that the solution of the 
exact Schrodingcr equation (|1.2|) with initial data in Bn,k for sonic finite K > N 
features an infinite rank at any positive time as long as many-body potentials are 
involved (see [TH] for related issues on the stationary solutions and Subsection 13.51 
for the picture for non- interacting electrons). We therefore have to rely on an 
approximation procedure that forces the solutions to stay on the set of ansatz for 
all time. In Physics' literature, the MCTDHF equations are usually (formally) 
derived from the so-called Dirac-Frenkel variational principle (see, among others, 
[T5| , IT6 l El] and the references therein) that demands that for all t £ [0, T], — 
£ B NtK and 



5*)=0, for all 5V £ T^dB 



N,Kt 



(3.1) (,•_-«* 

where T^Bn.k denotes the tangent space to the differentiable manifold 8Bn,k at 
iff. Equivalcntly, one solves 



(3.2) *(t) = argmiii|||i— -W*|| iB(0)TiL 2 ( n*r)) : * £ B N , K j 

for every T > (see [28]). A continuous flow t i— > *f?(i) £ &Bn,k on [0, T] may be 
lifted by infinitely many continuous flows 1 1— ► (C(t); <&(£)) foliating the fibers 8Tn,k 
that are related by the transitive action of a continuous family of unitary transforms. 
So called gauge transforms allow then to pass from one flow t i-> (C(t),$(t)) to 
another (equivalent) flow t h-> (C(t), $'(t)) such that #(i) = n(C(t), $(i)) = 
7r(C" (t), $'(£)). This is illustrated on Figure [fl below. 
One choice of gauge amounts to imposing 

(3.3) = for all 1 < i, j < K 

to the time-dependent orbitals. Formally the minimization problem ()3.2p under 
the constraints ^ = 7r(C, $), (C, $) £ Tn.k along with (|3.3| leads to the following 
system of coupled differential equations 

i r(c(t)) ^ = (i- p*) w^*[n *], 

(C7(0),$(0)) = (C ,$o), 
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for a given initial data (Co^o) i n ^ 7 n,k- This system will be referred to as the 
variational system in the following. 

The operator P$ in So denotes the projector onto the space spanned by the <^s. 
More precisely 

K 

(3.4) P*(-)=E<-,^>^- 
Actually one checks that 



Up to the Lagrange multipliers associated to (|3 . 3|) the right-hand side in the vari- 
ational system corresponds to the Frechet derivatives of the energy expectation 

£{^f) = (h \& | ^ with respect to the conjugate (independent) variables C and $. 

The variational system Sq is well-suited for checking energy conservation and 
constraints propagation over the flow as shown in Subsection 13. II below. However 
it is badly adapted for proving existence of solutions for the Cauchy problem or for 
designing numerical codes. Equivalent representations of the MCTDHF equations 
over different fibrations is made rigorous in Subscction l3.3l In particular, we prove 
below that the variational system is unitarily (or gauge-) equivalent to System (|3.26[) 
- named working equations - whose mathematical analysis in the physical case is 
the aim of Section [U 

a* 

Remark 3.1. Since for every a £ Sjy k, t; — = the system for the c a 's can 

dc a 

also be expressed as 

(3.5) i^- = J2(H$r\*<r)cr. 

T 

This equation is then obviously linear in the expansion coefficients. Furthermore, 
when the <f>i 's ( or equivalently the $ CT 's ) are kept constant in time, \3.5\) is noth- 
ing but a Galerkin approximation to the exact Schrddinger equation (|1.2p . The 
MCTDHF approximation then reveals as a generalization to a combination of time- 
dependent basis functions ( with extra degree of freedom in the basis functions ) of 
the Galerkin approximation. 

3.1. Conservation Laws. In this subsection, we assume the full-rank assumption 
on the time interval [0, T); that is Y(C(t)) is invertible for every t S [0, T]. Wc 
check here that the expected conservation laws (propagation of constraints, conser- 
vation of the energy) are granted by the variational system. Recall that to avoid 
technicalities all calculations in this section are formal but would be rigorous for 
regular classical solutions. We start with the following 

Lemma 3.2 (The dynamics preserves Tn,k)- Let (Co, 'I'o) € J 7 n,k being the initial 
data. If there exists a solution to the system Sq on [0,T] such thatr&nkT(C(t)) = 
K for all t G [0,T], then 

VMt)| 2 = l, / (f> i (t)fa{t)dx = 6 i , j , for allte[0,T}. 

Proof. First we prove that £ CT M*)| 2 = E CT M0)| 2 for all t e [0,T]. By taking 
the scalar product of the differential equation satisfied by C in Sq with C itself, we 
get 

j t \C{t)\ 2 = 2 Hj t C(t), C{t)) = 2 Q Y.( H * I = 2 Q(H * | *) = 0, 
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thanks to the self-adjointness of H, where 3? and 3 denote respectively real and 
imaginary parts of a complex number. From the other hand, the full-rank assump- 
tion allows to reformulate the second equation in (<So) as 

(3.6) = (i-p$)r(c)- 1 V***[W*]. 

(Notice that P$ commutes with F(C) -1 .) By definition I — P$ projects on the 
orthogonal subspace of Span{$}, therefore lives in Spanj^}- 1 for all t. Hence, 

(3.7) (*e> lWt) )_a 

for all 1 < i, j < K and for all t £ [0, T]. This achieves the proof of the lemma. □ 

We now check that solutions to the variational system indeed agree with the 
Dirac-Frenkel variational principle. 

Proposition 3.3 (Link with the Dirac-Frenkel variational principle). Let (C, $) G 
dJ~N,K be a classical solution to So on [0,T]. Then, \I r = n(C, <f>) satisfies the 
Dirac-Frenkel variational principle (|3.ip . 

Proof. We start with the characterization (|2.29|) of the elements in T^8Bn,k- Since 
the full-rank assumption is satisfied on [0, T], the orbitals satisfy (|3.6p . and therefore 
^ G Span!*}- 1 for all t G [0, T] and 1 < k < K. Firstly, being given a G E N , K , 
we have 

<"> (-f - Hf > - 'E t - (« +'S> <£l».> 

r r 

-«lM*«l».>-o. 

thanks to the equation satisfied by Co-. Indeed, 

dt ~ ^ d(j) k L at J 
fe= i 

and therefore the sum in (|3.8[) vanishes thanks to Lemma T2. 91 Secondly, for every 
1 < k < K and for any function ( in Span{$}- L , we have 

-<'(w£).-£™-0 

/ aw* \ 
(3.10) =-(p*— [w*],c)=o. 

Indeed, on the one hand, in virtue of Lemma l2.9l the first term in the right-hand side 
of (|3.9p vanishes whereas the second one identifies with % Ylfj=i ^kj(C) (^p, C) = 

((r(C(t)) ■ §|) , c) since ^ and C both belong to Span-}*}- 1 . On the other 

hand, the last line (|3.10p is obtained using the equation satisfied by <& in Sq and by 
observing that P$ ( = since ( G Spanj^}- 1 . The proof is complete. □ 

Let us now recall the definition of the energy 

£(*) = £(ir(C, $)) = (H*\*y 
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It is clear that depends on time via (C(t), <&(£))■ As a corollary to Proposi- 

tion [331 we have the following 



Corollary 3.4 (Energy is conserved by the flow). Let (C, $) € 8J-n,k be a solution 
to So on [0,T] such that 7r(C(t), $(f)) lies in the domain of Ti (or in the "form 
domain' when TL is semi-bounded) for all t in [0,T]. Then, 

£(n(C(t), $(t))) = £(n(C°, $ )) on [0, T]. 
Proof. Comparing with ()2.29p we observe that ^ e T^OBn.k, for 



K 



dt ~ M a + N 2^ 2^ C<T ftrh, 



cr k=l 



dt 



with in Span{$}- L whenever IT(i) is invertible. Then, applying Proposition 



dt 

to Sty = one obtains 



= sr - H * | ? \ = - St(wt I ?\ = -™0 



dt 



dt 



dt I 



2dt 



Hence the result. 



□ 



3.2. An a posteriori error estimate. We establish an error bound in L 2 (Q) N 
for the MCTDHF approximation compared with the exact solution to the linear 
TDSE (|1.2p . Let us introduce the projection Vt^oBn^k on t° the tangent space 
T^8Bn,k to 8Bn,k at \S r . Then, we claim 

Lemma 3.5. Given an initial data (C°, <fr°) € 8Tn,k and an exact solution \& e to 
the N -particle Schrodinger equation (|1.2|) . Then, as long as (C, <&) zs a solution to 
S in dFx,K, we have for = 7r(C(£), $(£)) and \I> = 7r(C°, $°) £/ie estimate: 

||^-*|| i2(nJV) < || 1 M0)-*°|| i 2 (n * ) +y ||(J-7*r.a8*,J[«*(s)]|| La(n) * ^. 

Proof. First, Proposition 13.31 expresses the fact that Vt^oBn k (*'lif — = 0. 

Therefore the equation satisfied by the ansatz "J is: 



9i 



(3.12) i—-H^ = (I-P T ^B N , K ) i^r-n* =-(I-V T *8B N , K )[H*] 



dt 



since ^ lives in the tangent space T^8Bn,k- Next, subtracting (13.12p from ([T 
we get 

(3.13) ^ E m ^ *) = -(I VtoBB*,*) [H *] 

Then, we apply the PDE above to ^>e — and we integrate formally over fl N . The 
result follows by taking the imaginary of both sides and by using the self-adjointncss 
of H. □ 



Roughly speaking, the above lemma tells that the closer is Ji "J to the tangent 
space T^OBn.k, the better is the MCTDHF approximation. Intuitively, this is true 
for large values of K. Let us mention that this bound was already obtained in [28] 
and it is probably far from being accurate. However if the MCTDHF algorithm is 
applied to a discrete model say of dimension L then for K large enough (K > L ) 
this algorithm coincides with the original problem (see Subsection 17. 3p . 
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3.3. Unitary Group Action on the Flow. The variational system So is taylor- 
made for checking energy conservation and constraints propagation over the flow. 
However it is badly adapted for proving existence of solutions for the Cauchy prob- 
lem or for designing numerical codes. It is therefore convenient to have at our 
disposal several explicit and equivalent representations of the MCTDHF equations 
over different foliations and to understand how they are related. This is the purpose 
of this subsection. Proofs of technical lemma and theorems are postponed in the 
Appendix to facilitate straight reading. 

We start with the following (straightforward) lemma on regular flows of unitary 
transforms : 

Lemma 3.6 (Flow of unitary transforms). Let Uq £ Ok and let t i— > U(t) be in 
C x {[Q,T);O k ) with U(0) = U Q . Then, 1 1-> M(t) :=-i^-U defines a continuous 
family of K x K hermitian matrices, and for all t > 0, U(t) is the unique solution 
to the Cauchy problem 

, ^ \i^T~ = U(t)M(t), 

(3.14) I dt w w ' 

{ U(0) = U . 

Conversely, ift*—> M(t) is a continuous family of K x K Hermitian matrices and 
if Uq £ Ok is given, then {3.1$ defines a unique C 1 family of K x K unitary 
matrices. 

The corresponding flow for unitary transforms on expansion coefficients is as 
follows: 

Corollary 3.7. Let (N,K) be an admissible pair, let t t— > M(t) be a continuous 
family of K x K Hermitian matrices and let U° £ Ok ■ Then, ift^U (t) denotes 
the unique family of unitary K x K matrices that solves (|3.14[) . the unitary r x r 
matrix U given by (|2.16p is the unique solution to the differential equation 

(3.15) ^ dt 



\ U(0) = d(U°) 



with 

(3.16) M^ T = V (-lJ^W-H-^Af, 



/ 4 V ''I "^3- 
<A«=t\{j} 



The proof of this corollary is postponed to the Appendix. The main result of 
this section is : 

Theorem 3.8 (Flow of unitary equivalent foliations). Let Uq £ Ok and (Co, $o) € 

(i) Let t i — ► M(t) be a continuous family of K x K Hermitian matrices on [0, T] 
and let U(t) £ C 1 ([0, T); Ok ) be the corresponding solution to (|3.14p . Assume that 
there exists a solution (C, $) £ C°(0, T; 8Tn,k) °f So with initial data (Co,3?o)- 
Then, the couple (C, $') = U{t) ■ (C, $) with U £ O r K defined by ([2H5J) and 
is solution to the system 

( dC I \ 

*^- = (w*|Vc*) -M'C, 

(3 - 17) 1 iT(C) ^- = (I-P*0 V*/**[H*]+r(C)M / * / 

dt 

t (C(0),$'(0)) = %-(C ,$o) 
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with f = tt(C*,$) = n(C',<f>'), U = (U ,d(U )) G O r K being defined by and 

M' = UMU*, M' = UMU*, 

where M is the r x r Hermitian matrix with entries given by (|3.16|) . 
('iij Conversely, assume that there exists a solution (C, 3?) G C° (0, T; dT^ t K ) to 
5o wi£/i initial data (Co,$o) arl ^ fef ^(^) G C 1 ([0, T 1 ); 0^ ) . Then, the couple 
(C',&) = U(t) ■ (C,$) with U G 0^ denned 6?/ (pTT%|) and is a soZuiion to 

%stom (j3~T7|) with M (t) = -i U . 



Remark 3.9 (Link with Lagrangian interpretation). The equations can be derived 
(at least formally) thanks to the Lagrangian formulation: One writes the stationar- 
ity condition for the action 



Am = f 

Jo 



i— WkVk)dt 

at 1 



over functions iff = iff (t) that move on J^n,k- The associated time- dependent Euler- 
Lagrange equations take the form (|3.17p with iff = 7r(C, $), M an hermitian matrix 
and with M be the rxr hermitian matrix linked to M through Eqn. (|3.16[) above. As 
observed already by Cances and Le Bris [5], even if they appear so, the Hermitian 
matrices M and M should not be interpreted as time- dependent Lagrange multipliers 
associated to the constraints (C, $) G J-n.k since the constraints on the coefficients 
and the orbitals are automatically propagated by the dynamics (see Lemma \3.15\) . 
but rather as degrees of freedom within the fiber at iff . Ln particular, this gauge 
invariance can be used to set M and M to zero for all t as observed in Lemma \3.8\ 
and Eqn. (|3.23p below, so that the above system can be transformed into the simpler 
system (Sq ) we started from. 

As a first example of the change of gauge one can use the unitary transforms to 
diagonalizc the matrix F for all time and therefore derive the evolution equations 
for natural orbitals following [5] 

Lemma 3.10 (Diagonal density matrix). Let (C, $) satisfying Sq with initial data 
(Co,$o) an d let Uq G Ok that diagonalizes F(Co). We assume that for all time 
the eigenvalues ofF(C) are simple, that is 7, 7^ jj for 1 < i,j < K and i 7^ j. 
Define a K x K Hermitian matrix by 



1 



Ma, 



7 j - It 




dif 



d<f>j 



if i + 3, 
otherwise 



and consider the family t t— > U(t) G Ok that satisfies (|3.14p with U(t = 0) = Uq. 
Then (C, $') = U{t) ■ (C, $) is solution to 



dC 
~dT 



1 7 ,(t) M = (I _ Pft ,) ^L\n *] + 7l .( f ) M' $' 

^ (C"(0),$'(0)) = W -(Co,$o) 

wii/i i/ie notation of Theorem \3.8[ Ln particular, F(C) = diagyy\(t), . . . ,7r- (t)J 
/or every t. 
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Proof. Using the equation for the coefficients in (|3.17[) together with (|2.1ip . the 
evolution equation for the coefficients of the density matrix writes 

, drf i: 



1?= E <-•>■ 



^ [{H # | c a * T ) - (c r $ CT | ft *) 



cr,r : j'Gt 
<AW=t\{j} 

E 

<A{i}=r\{j} 



(-l) c 



ft* 



ft* 



if 

k=l 



Mifc r 



Next, we require that 

7i?M = 7* that is 



V 1 < i ^ j < K. 



Using the above equation, a sufficient condition is given by 



Mi, 



1 



H ~ 1] 



ft * 



ft * 



□ 



This achieves the proof. 

As a second application of Theorem 13.81 we investigate particular (stationary) 
solutions or standing waves. A standing wave for the exact Schrodingcr equation is 
of the form &(t, x) = e~~ lX t ^(x) with A 6 R . In the same spirit we look for solutions 
(C, *') of System (l3~T71) with (C, $') = U{t) ■ (e~ iXt C, $), where (C, $) G dF N ,K 
is fixed, independent of time, and U(t) 6 O^-. Using the formulas (|2.25|) and (I2.26P 
for the changes of variables, we arrive at 

/. dU(t) 



dt 



MU + AU C 



T(C) (i U* 



dU 
~dt 



P*) V*** [ft *], 



U* M U) $ = (I - 
17(0) = Ik-. 

In the above system \I/ = 7r(C, $) and T(C) are independent of time and T(C) is 
invertible. We start with the equation satisfied by <f>. Observing that the left-hand 
side lives in Span{$} whereas the right-hand side lives in Span{$}- L , we conclude 
that there are both equal to zero. Therefore, there exists a K x K matrix A that 
is independent of t and such that 

V$**[ft *] = V ¥ (ft * | *\ = A • $. 



(3.18) 

Also since the left-hand side has to be independent of t we get 

dU 

i — = MU. 
dt 

Comparing now with the equation for the coefficients we infer from Corollary 13.71 
that 

. dll{t) 



hence 
(3.19) 



dt 



( ft * I * 



-MU, 



AC. 



Equations (|3 . 18[) and (|3 . 19[) are precisely the MCHF equations that are satisfied by 
critical points of the energy. They were derived by Lewin |26j in the Coulomb case. 
The real A is the Lagrange multiplier corresponding to the constraint C <G S r ~ 1 
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whereas the Hermitian matrix A is the matrix of Lagrange multipliers corresponding 
to the orthonormality constraints on the orbitals. Existence of such solutions in 
physical case is recalled in Section |U 

The proof of Theorem 13.131 is postponed in the Appendix and we rather state 
before some corollaries or remarks. In Physics' literature the MCTDHF equations 
are derived from the variational principle (|3.1|) under the constraints <I> = ir(C, G 
Bn.k along with additional constraints on the time-dependent orbitals 

(3.20) <f>j) = <G0 4 , <f>j) for all 1 < i, j < K. 

In the above equation G is an arbitrary self-adjoint operator on L 2 (fl) possibly 
time-dependent named the gauge. In this spirit the variational system corresponds 
to G = 0. Therefore a gauge field is chosen a priori and the corresponding equations 
are derived accordingly. Both approaches are equivalent by observing that, to every 
Hermitian matrix M, one can associate a self-adjoint operator G in L 2 (Q) such that 
Mij = {G<f>i,<pj) by demanding that 

K 

G 4>i = ^2 M ij <t>3 for au 1 < i < K - 

3=1 

Conversely being given the family t i— > M{t) in Theorem l3.13l it follows immediately 
from the system (|3.17[) that for all 1 < i,j < K, 




provided T(C) = UT(C) U* is invertible on [0,T). We state below Theorem OH] 
that is the equivalent formulation of Theorem 13. 131 in terms of gauge. It is based 
on above remarks together with the following : 

Lemma 3.11. Let t i— > G(t) be a family of self-adjoint operators on L 2 (fl) and let 
<t> = (<f>i(t), <f>2(t), . . . , 4>k if)) € Cl 2 (o) k such that such that t i-> (G(t)^(t), <j)j{t)) 
is continuous on [0,T] . Then the matrix M with entries Mij(t) = (G(t)4>i(t),4>j(t)) 
is Hermitian and the Cauchy problem (|3.14[) defines a globally well-defined C 1 flow 
on the set of unitary K x K matrices. In that case, the unitary transforms U = d(U) 
solve the Cauchy problem (|3 . 1 5[) with M in (|3.16p given by 

N 

(3.21) M^ T =Y / (G X ^a\^r)- 

i=\ 

Remark 3.12. In Lemma \3.11\ the functions t t— > 4>i(t) are continuous with val- 
ues in the domain of G . When G is bounded from below it is enough to assume 
continuity in the form- domain. When G is the Laplace operator or, more gen- 
erally a one-body time-independent Schrddinger operator, we simply assume that 
4>,i G -ff^R 3 ) or 4>i G Hq(VI) when Q is a bounded domain. (Other boundary condi- 
tions could of course be considered.) 

Theorem 3.13 (Flow in different gauge). Let Uq G Ok, (Co^o) G 8J-n,k and 
let t h G(() k a family of self-adjoint operators in L 2 (fl). Assume that there 
exists a solution (C, $) G C°(0, T; OFn.k) to So with initial data (Cq,^>q) such 
that t h- > (G(i)0i(t), 4>j(t)) is continuous on [0,T] for every 1 < i, j < K. Define 
the family of unitary transforms U{t) G C 1 ([0, T); Ok ) that satisfy (|3.14[) with 
Mi 3 = {G4>i,4>j) as in Lemma [JU[ Then the couple (C",<I>') =U(t) ■ (C, $) with 
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U{t) = (d(U(t));U(t)) defined by (f2~T5|) and (gTTgD is a solution to 

N 



iT{C'\ 



dt 



N 



= F(C') G & + (I - P#,) V*-**[W * - ^ Gs, . 



I (C(0),$'(0)) = Wb-(C ,*o), 

wit/i * = tt(C*,$) = tt(C",$')> = (U ,d(U )) G 0^ H denned 6j/ (l2~18l) and 
iw'i/i M being the r x r Hermitian matrix given by (|3.16[) . 

Remark 3.14. Passing from So to Sg amounts to change the operator H by H. — 
Y^Li G Xi in both equations and by adding the linear term F(C")G$' in the equation 
satisfied by Note that whereas solutions to So in OJ-n.k satisfy 

for all 1 < A', solutions to (Sg) satisfy 

(3-22) i<^>«> = <GM>. 

System Sg corresponds to the choice of gauge G . 

This is illustrated and and summarized on Figure [T] below. 

Theorem 13. 131 and Lemma 13.111 provide with the differential equation that satis- 
fies the unitary matrix U(t) that transforms Sq into Sg- A direct calculation shows 
that, given two self-adjoint one-particle operators G and G, the solution to 



(3.23) 



■dU 

i-jr = U M G ^ G , 
dt 



U(t = 0) = U° 

with (Mc^g 1 ) = ((G — G')4>i, 4>j) maps a solution to Sg to a solution to Sg' ■ In 

particular, if we prove existence of solutions for the system Sg for some operator G 
then we have existence of solutions for any system 5g' • Another immediate though 




Figure 1 . Flow on the Fiber Bundle 
crucial consequence of Theorem 13.131 and Theorem 13.81 is given in Corollary 13.151 
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below. It states that for any choice of gauge the constraints on the expansion 
coefficients and on the orbitals are propagated by the flow and the energy is kept 
constant since it is the case for the system So- Also the rank of the first-order 
density matrices does not depend on the gauge. 

Corollary 3.15 (Gauge transforms and conservation properties). LetT > 0. LetG 
be a self-adjoint (possibly time- dependent) operator acting on L 2 (fl). Assume that 
there exists a solution to the system So on [0, T] such that rankT(C(<)) = K and 
such that the matrix t t-y (G(j>i, <^j) 1<i j <K is continuous. Then, for all t G [0,T], 

(C(i), $(t)) G dF N , K , 
and the energy is conserved, that is 

In addition, $ = 7r(C, <£>) satisfies the Dirac-Frenkel variational principle (|3.ip . 

Proof of Corollaru VS .15\ By Theorem l3.13l and its remark, if (C, $) satisfies Sg with 
initial data in 8Tm,k, there exists a family of unitary transforms U G C 1 (0, T; Ok) 
such that (C, $) = U ■ (C",$') where (C",$') satisfies Sq with same initial data. 
Since by Lemma l3~2"l So preserves Tn,k, so does iSg since U and U = d(U) are 
unitary. Then, by Lemma l2~5l 7r(C, $) = 7r(C", $') = ^, and the energy is conserved 
by the flow since it only depends on \E'. Eventually Eqn. (|3.1[) is satisfied since 
T^dBjq.K only depends on the point ^ on the basis Bn,k and not on the pre- 
imagcs in the fiber 7r _1 (^l/). □ 

So far we have considered a generic Hamiltonian TL and we have written down an 
abstract coupled system of evolution equations for this operator. In the following 
subsection we turn to the particular physical case of TV-body Schrodingcr-type 
operators with pairwise interactions 



3.4. iV-body Schrodinger type operators with pairwise interactions. At 

this point, we consider an Hamiltonian in L 2 (VL N ) of the following form 



iV 



(3.24) 



i=l l<i<j<N 



In the above definition, H is a self-adjoint operator acting on L 2 (fl). To fix ideas 
we take H = — + U. v is a real- valued potential, and we denote 

V= v(\xi-Xj\). 

l<i<j<N 

Expanding the expression of TL in the system So and arguing as in the proof of 
Theorem 13. 131 we obtain 



(3.25) 



dC 



N 



So ■ < 



dt 
<9$ 



N 



iT(C) — = (I-P*) V^*[V^ + Y, H Xj 
I (C(0),$(0)) = (C >*o) &Fn,k- 
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Comparing with System Sg in Theorem l3.131 one observes that the choice of gauge 
G = H leads to the equivalent system 

. dC 



(3.26) 



5; 



H 



i T{c) ^ = r(c) h $ + (i - p») v$**[y *] 



^ (C(0),$(0)) = (C ,$o) £ J~. N.K: 

(provided t 1— > (H.<j>i,(j>j) makes sense). From Corollary 13.151 we know that if the 
initial data in (|3.26p lies in J~n.k it persists for all time. This property allows to 
recast System (|3.26j) in a more tractable way where the equations satisfied by the 
orbitals form a coupled system of non-linear Schrodinger-type equations. This new 
system that it is equivalent to System (|3.26|) as long as the solution lies in J~n,k 
will be referred to as working equations following [JJ 124| . It is better adapted for 
well-posedness analysis as will be seen in the forthcoming section. 

Proposition 3.16 (Working equations). Let (C, $) be a solution to (|3.26p inTN,K, 
then it is a solution to 



(3.27) 



i r(c) ^ = r(c) h $ + (i - p«) w[c, $] $, 



[ (C(0),$(0)) = (C°, $°) G J~N,K 1 
where K[3>] (resp. W[C, $]J is a r x r (resp. K x K) Hermitian matrix with entries 
(3.28) 



and 
(3.29) 



A' 



W[C,*] <J .( a? ) = 2 2 
fe,;=i 

where here and below we denote 



D v (f,g) 



v(\x-y\)f(x)g(y)dxdy, 



nxn 



?'(• -y)f(y)dy 



and with the coefficients Jijkl being defined by (|2.10[) m Provosition \2.3\ Conversely, 
any solution to (|3.27p defines a flow on Tn,k as long as T(C) is invertible and is 
therefore a solution to (|3.26|) . 

Proof. We have to show that for \I> = tt(C, $) in S/v^a 

(3.30) (V" * I V c $) = (y # I *) = K[i>] C 
and 

(3.31) V^*[V = v ¥ (V * 1 *> = W[C, $] 
We start from 



(y * 1 *) 



«) *] :2 (x, y, x, y) - y|) <&rofr/ 



with 



A 



^^} :2 (x,y,x,y) = T^fei <t>i{x) <t>j(y) <t>k( x ) Mv) 
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according to (|2.9[) . Since only the coefficients jijki depend on C through Eqn. 
(gUl we first get 

K 

\7c(V * | f ) = V c (7tffcl) A, ^ . ■ 

i,j,k,l— 1 

Hence (|3.28|) by using again Formula l|2.10[) . 

We now turn to the proof of (|3.31[) starting from 

(Vtf|tt) = V Jijki ff (j) l (x)(t) J {y)^ k (x)^ l (y)v{\x-y\)dxdy. 

Then, for every 1 < p < K 

8 K - K 

-=- (V | *) = ^] 7 ypZ ((fa<h)*v) fa + J] 7 4 jfc P ((0<0*)*«) 0j 

^P i,j,i=l i,i,fc=l 

if 

by interchanging the role played by i and j in the first sum and by using jijkp = "fjipk 
and renaming fc as / in the second one. Comparing with (|3.29p we find 

8 K 

(V* |*) = 2^W[C,$] B -^. 

9 ^p j=J 

To achieve the proof of the proposition we now check that the system of equations 
in (|3.27p preserves Tm,k a s long as T(C) is invertible. The claim is obvious as 
regards the orthonormality of the orbitals since H is sclf-adjoint and since I — P$ 
projects on Spanj^}- 1 . On the other hand, the equation on the coefficients leads 
to 

^||C(t)|| 2 = 2 3^K[$] CT . rCrC - CT = 

(T.T 

since the matrix K[$] is Hermitian. □ 



We treat apart in the last two subsections the special cases of the linear free sys- 
tem with no pairwise interaction and of the time- dependent Hartree-Fock equations 
for the evolution of a single-determinant (TDHF in short) with pairwise interaction. 

3.5. Interactionless Systems v = 0. In this section we consider systems for 
which the binary interaction potential v is switched off. Then the system (|3.26l) 
becomes 

.dC 

l H = °' 

iT(C) ^ =T(C)H$. 

From the first equation the coefficients c^'s are constant during the evolution. In 
particular the full-rank assumption is satisfied for all time whenever it is satisfied 
at start. In the latter case the orbitals satisfy K independent linear Schrodinger 
equations 

(3.32) Z _ =H $, 
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and the -ZV-particle wave-function $ = tt(C, $) solves the exact Schrodinger equa- 
tion 

*(f = 0) =^(C ,$o). 

Conversely, the unique solution to the Cauchy problem (|3.33[) with (Co,$o) G 
dJ^N,K coincides with 7r(C(f), <&(£)) € Fn.k where is the solution to (|3.32j) . 
This is a direct consequence of the fact that the linear structure of (|3.33|) propagates 
the factorization of a Slater determinant. In particular, this enlightens the fact 
that the propagation of the full-rank assumption is intricately related to the non- 
linearities created by the interaction potential v between particles. 

3.6. MCTDHF (K = N) contains TDHF. The TDHF equations write (up to 
a unitary transform) 

(3.34) i^=H4>i+F 9 <t>i, 

for 1 < i < N, with J 7 ® being the self-adjoint operator on L 2 (fl) that is defined by 

N N 

F$w=(^2 v {\- -y\)\<t>3{y)\ 2d v) w ~^2\ / v (\ ■ -v\3j{v)w{y) dyj fa. 

j — 1 j — 1 

The global-in-time existence of solutions in the energy space 77 1 (J7) Ar goes back to 
Bove, Da Prato and Fano [B] for bounded interaction potentials and to Chadam 
and Glassey [T^] for the Coulomb potentials. They also checked by integrating the 
equations that the TDHF equations propagate the orthonormality of the orbitals 
and that the Hartree-Fock energy is preserved by the flow. Derivation of the 
TDHF equations from the Dirac-Frenkel variational principle may be encountered 
in standard Physics textbooks (see e.g. [3D]). Let us also mention the work [5] by 
Cances and Le Bris who have investigated existence of solutions to TDHF equations 
including time-dependent electric field and that are coupled with nuclear dynamics. 

By simply setting K = N in the MCTDHF formalism one gets 

(3.35) = 1 , F(t) = I N 
and 

:= C{t) 0i (i) A ... A c/> N (t), C(t) = e" ie *W 
for some 0$ € R. In addition according to Remark 12.41 



(3.36) jjka = - {pi,j 5 k ,i - 6i t k 5 3 ,i)- 

Therefore with the definitions (|3~2"51) and ([3~2^]) 

K[$] = ( V(f>i A . . . A (j) N | fa A . . . A 4> N ) 

i,j,k,l : {i,j} = {kd} 
N 

= 2J (F$fa,fa) 



and 



N 

X)W[C,$]y ^ = T*fa 
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Eventually for N = K, according to (|3~?7|) . the MCTDHF system in the working; 
form turns out to be 



N 



S n (N = K) { 



dO$(t) 
dt 

i^ = H0i + (I-P*)^«, 

JV 



= H 4>% + & - fyi^j) 4>j 

with #$(0) = and $(0) 6 O^^n . Comparing with (|3.17|) . we introduce the N : 
N Hcrmitian matrix M with entries My = — ( ^$ 0^ , <fij ) . According to Lemma 
there exists a unique unitary matrix U(t) such that 

dt 

U(t = 0) = Eat. 



In virtue of (|2. 16|) the unitary matrix U that transforms <fii A . . . A 4>m into (Utfii) A 
. . . A (U4>n) is then simply a complex number of modulus 1 (U = det(t/)) that 
satisfies 



(3.37) 



dU 

U(t = 0) = 1 



tr(M) U, 



Comparing (|3.37[) with the equation satisfied by 9<s>(t) in Sn{N = K) it turns out 
In that special case a change of gauge is simply a multiplication 



that U 



by a global phase factor. Applying Theorem 13.81 the functions 1 < i < N, 
defined by <&' = satisfy the standard Hartree-Fock equations (|3.34p and C'(t) = 
UC(t) = 1 for all time; that is ^ = <fr[ A ... A 4>' N . Being a special case of the 
MCTDHF setting we then recover "for free" that the TDHF equations propagate 
the orthonormality of the initial data, that they satisfy the Dirac-Frcnkel variational 
principle and that the flow keeps the energy constant. 



4. Mathematical analysis of the MCDTHF Cauchy Problem 



This section is devoted to the mathematical analysis of the Cauchy problem for 
the A-body Schrodinger operator with "physical interactions" 



(4.1) U(x) 
that is given by (|3.2GD : 



E 



x — R„ 



v{x) = R ' 



Si 



H 



dC I 

i t(c) — = r(c) h $ + (I - p#) v***[y 



k (C(0),*(0)) = (C 0) *o) eT N 



K ■ 
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In this section £1 = R 3 . According to Proposition (|3.27p solutions to (6>h) lie in 
Fn,k and they are therefore solutions to 

dC 



(4.2) 



with 



i T(C) = TOO) H $ + (I - P«) W[C, *] 

[ (C(0),$(0)) = (C°,$°) GJfjv.K 



K i 

W[C,$] ij (x) = 2 HUl{<l>k4>l*T-7) 
k,t=l ' ' 

D(f, g) = [[ — ^ f(x) g(y) dxdy . 

J Jr 3 xr 3 \ x _ y\ 

The above system is referred to as the "strong form" of the working equations. Let 
us emphasize again that it is equivalent to (<Sh) provided (C, $) G Fn,k- The main 
sources of difficulties arise from the fact that the matrix T(C) may degenerate and 
from the Coulomb singularities of the interaction potentials. Our strategy of proof 
works for more general potentials U and v. This is discussed in Section [7] below. 

The spaces C r and W m ' p (E. 3 ) K are equipped with the Euclidian norms for the 
vectors C and $ respectively 



K 



||C|| 2 := J2 \^ II*IIU-:=E IL 
Moreover, for a p x p matrix M we use the Frobenius norm 



lAsrii = 



E \ M > 



2 



We introduce the spaces X m := C x H m (R 3 ) K for to G N endowed with the norms 

ll(C»*)llx m = ||C|| + ||*||H-.. 
The main result in this section is the following 

Theorem 4.1. [The MCTDHF equations are well-posed] Letm>\ and (C°, $°) e 
d^N,K with <£>° m .//"(R 3 )^. Then, there exists a maximal existence time T* > 
(possibly +oo 6ui independent of m) such that: 

(i) The MCTDHF system (|3.27p admits a unique solution (C, $) luii/i 

C G C 1 ([0, T*); C r ) , $ G C°([0, T*); ff m (R 3 )*) n C* 1 ([0, T*); if m_2 (M 3 ) K ) . 

This solution depends continuously on the initial data (C , $ ) in X m . For every 
0<t<T*, 

(ii) (C(t),$(t)) G dF NtK and T(C(t)) is invertible. 
(Hi) The energy is conserved : 

(ttjv *(*)|#(t)) = (WArf !* ) Mi/i # = 7r(C,$) and *° = 7r(C°,$°). 

(iv) The Dirac-Frenkel variational principle (|3.ip zs satisfied. 
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(v) When T* < +00 one has 



limsup||r(C(t))- 

tyr* 



and more precisely: 



T* 



n 



\\T(C(t)) 1 \\ 3 / 2 dt = +oo. 



The global well-posedness in H 1 and H 2 of the TDHF equations goes back 
to Chadam and Glassey [12]. Recently Koch and Lubich [24] proved local well- 
posedness in H 2 of the MCTDH and MCTDHF equations for regular pairwise 
interaction potential v and with U = by using Lie commutators techniques. Our 
result extends both works. The rest of the section is devoted to the proof of this 
theorem. The above system with the same notation is rewritten in the "mild form" 
which makes sense as long as the matrix F(C(t)) is not singular: 

(4.3) U(t) = e- itA U - i f e- t[t - s)A C(U{s)) ds 

Jo 

with 
(4.4) 

_ / C \ / \ rnn _f K[*]C 

V $ ) \ H®Ijc ) ' [ ' ~ \ T(C)- 1 (I-P*)W[C,$]$ 

The strategy of proof is as follows. 

In Subsection 14.11 we show that the operator L is locally Lipschitz continuous 
on X m for to > 1 in the neighborhood of any point (Co,$o) such that F(Co) is 
invertible. Observe in particular that T(C) is a second-order homogeneous function 
of the coefficients C and therefore the invertibility of this matrix is a local property. 
Standard theory of evolution equations with locally Lipschitz non-linearities then 
guarantees local-in-timc existence and uniqueness of a mild solution in these spaces 
that is continuous with respect to the initial data as long as the matrix F(C) 
remains invertible (see e.g [331 |32 ). Next for initial data in X m with to > 2. 
the corresponding mild solution in this space is regular enough to be a strong 
solution to (|4.2p (see [32j[TT]). As shown in the previous section (Proposition [3TG|, 
the strong solution then remains on the constraints fiber bundle OJ-n.k and it 
is therefore a solution to (|3.26[) . Furthermore using the gauge equivalence and 
Corollary 13.151 one deduces that the energy of the solution is conserved and that 
the Dirac-Frenkel variational principle is satisfied. Recall for further use that the 
energy may be recasted in the following equivalent forms |17[ 126] 

(4.5) 

£(tt) = £(ir(C,$)) 



h r + -w[c, $] 



L 2 (o.y 



K 



iv<^ • + u <j)i <f>j 



K 

dx + 2J lijki D (<pi 4>i ; fa (f> 3 1 

ZJjfc,/ — 1 



In consequence for initial data in X m , to > 2, the norm of the vector <E>(i) remains 
locally bounded in H 1 (independently of the H 2 norm). Therefore it is also a 
strong solution in H 1 defined on the same time interval which depends only on the 
i? 1 norm and on T(Co). Eventually using the density of X2P\dJ 7 N,i< in X\ nd!pN,K 
and the continuous dependence with the initial data one obtains the local-in-time 
existence of a strong solution in X\ D 8Tn,k with constant energy. 
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In Subsection 14.21 relying on the conservation of the energy we prove the ex- 
istence of the solution over a maximal time interval beyond which the density 
matrix degenerates. The equations themselves imply the further regularity C(t) £ 
C 1 ([0,T'),C7) and$(t) e xC°([0, T*), H m (M. 3 ) K ) n C 1 ([0, T*), H m ~ 2 (M. 3 ) K ) . 

4.1. Properties of the one-parameter group and local Lipschitz properties 
of the non-linearities. As in Chadam and Glassey [T^] for example, one checks 
that {e ltA } is a one-parameter group of linear operators, unitary in Xq and 
uniformly bounded in time for < t < T in X\ and Xi. 



We now show that the operator C in the right-hand side of (J473J) is a locally 
bounded and locally Lipschitz continuous mapping in a small enough neighbor- 
hood of any (Co,^»o) in X m such that F(Co) is invertible for every to > 1. The 
operator C reveals as a composition of locally bounded and locally Lipschitz con- 
tinuous mappings as now detailed. We first recall that invertible matrices form an 
open subset of A4kxk(C) and that the mapping M i— > M^ 1 is locally Lipschitz 
continuous since 

p/ _1 - M~ 1 \\ = WM- 1 (M - M) M" 1 !! 

< ||M _1 || ||M _1 || ||M-M||. 

In addition, being quadratic, the mapping C > F(C) is for any to and indepen- 
dently of to locally Lipschitz in X m in a small enough neighborhood of any (Co, $0) 
such that T(Co) is invertible. The same holds true for the mapping C 1— > T(C) -1 
by composition of locally bounded and locally Lipschitz functions. 

The operator P$ is a sum of K terms of the form ((/>, -)ui <fi with <f> in iJ Tn (R 3 ). 
Hence, for m > 0, 

(4-6) I|P*IU(1P») Z PIU 2 

where here and below < is a shorthand for a bound with a universal positive 
constant that only depends on K and N. Therefore $ 1— ► P$ is locally Lipschitz 
from X m to C{H m ) since it is quadratic with respect to <E>. To deal with the other 
non-linearities we start with recalling a few properties of the Coulomb potential 
taken from [121 Lemma 2.3]. Their proof is a straightforward application of Cauchy- 
Schwarz' and Hardy's inequalities and we skip it. Let </>, tp <E i? 1 (R 3 ), then with 
r = \x\, (0V) W 1 ' 00 ^ 3 ), and we have 

(4.7) ll(#)*;IU-(R»)<2||V0|| £ 2 (n) M| £a( n) 
and 

V((^)*i) <4||V0|| ia( n)||V^|| L3( n). 

As a consequence of above inequalities and by an induction argument that is detailed 
in [H] for example, we have, for $ £ H m (M. 3 ) K and for every 1 < i, j, k < K, 

(4-8) ||((&&)*;)&L»CR3) £ !l*ll 2 Hm ' < ll^ll 3 ^ 

with to' = max(m— 1,1). First, recall from (|3.29p . that W[C, <&] $ is a sum of terms 
of the form •fjku 4>j (^*4>k </>z) with the coefficients jjku depending quadratically on 
C according to (|2.10[) . They are therefore locally Lipschitz continuous with respect 
to C. Gathering with (|4.8p we have 

(4.9) ||W[C,$] <f\\ Hm < ||C|| 2 ||$||^ m , \\$\\ Hm < \\C\\ 2 \\<£\\ 3 Hm . 

The mapping (C, $) 1— > W[C, $] $ is then locally bounded in X m and being qua- 
dratic in C and cubic in <E> it is locally Lipschitz continuous in X m by a standard 
polarization argument. In particular, the first bounds reveals a linear dependence 
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on the H m norm. Eventually, for every 1 < k, I < K, using (|4"7T)) and Holder's 
inequality we obtain 

\d(^^ i: M t )\ < 11*111, n$ii fl i <\m\% m , 

the last line being a direct consequence of (|4.8j) . In particular this proves 
(4-10) |K[*]|< ffli a ||*||Hi, 

(4.H) l|K[*]C||<||C||||$|| 4 ffm 

and that (C, $) i— ► K[*] C is locally Lipschitz continuous in X m since according to 
(|3.28p . K[$] C is a finite sum of terms of this kind up to some universal constant. 

For any m > 1 existence and uniqueness of a solution (C(t), *(i)) to the integral 
equation (|4.3p in a neighborhood of (C°, $°) in C°(0, T; X m ) for < T small enough 
follows by Segal's Theorem |33) . which also ensures the continuity with respect to 
the initial data in X m . 

We now turn to the existence of a maximal solution and to the blow-up alterna- 
tive in X\. 

4.2. Existence of the maximal solution and blow-up alternative. To sim- 
plify notation, from now on we use the shorthand T(t) for T(C{t)). Existence of 
a global-in-time solution requires to control uniformly both the H 1 norm of <& and 
the norm of T -1 (t). With the conservation of the energy this turns to be equivalent 
to control only the norm of T -1 (t). Let T* denotes the maximal existence time 
and assume that T* < +oo. Wc first show that 

(4.12) limsuplirfXT 1 !! = +oo. 

We argue by contradiction and assume that there exists a positive constant Mq 
such that for all t e [0,T*), ||r(t) _1 || < Mq. We now prove that there exists a 
positive constant Kq such that 

(4.13) Vie [0,7*), ||*(i)|| H i <K . 

Thanks to Lemma 13.41 and Corollary 13. 151 the energy is preserved by the flow, and 
therefore using the expression (|4.5[) 

(h r(t) *(t), *(*)) < (h r(t) *(<), $(*)) + \ (w[c, *] *(t), *(*)) 

= £(tt(C,*)) = £(.(C°,$ )) 
for all < * < T* since, with * = tt(C, *), 

{Vl-j^ = (W[C,*] $(*),$(<)) > 

for w > 0. As in [5H1[T7], Kato's inequality then yields that 

l|Vr*||Hi < Mi 

where Mi is a positive constant independent of t > 0. Now let <G (0, 1] be the 
smallest eigenvalue of the hermitian matrix IT(f). Then according to the definition 
of the Frobenius norm 

lir- 1 ^)!! = 



with {/ii(t), • • • , nxit)} being the eigenvalues of F(t), hence 



< V < ^tt and , < VT" 1 < 
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for all te[0,T*). Therefore 

(4.14) \\<£\\ m < l|VT$|| H1 < KV'Mx WT-'W 1 ' 2 . 

In particular, this shows (|4. 13|) with Kq = K l / A MiM^ 2 . Therefore, for any 
t G [0, T*) arguing as above, we may build a solution to the system on [t, t + T ] for 
To > that only depends on Mq and Kq. Since t is arbitrary close to T* we reach 
a contradiction with the definition of T* . Hence (|4.12l) . 

Now, taking the derivative with respect to t of both sides of IT F _1 = Ik, we get 

(4.15) — — = -r- 1 — r _1 , 

y ' dt dt 

for all t G [0,T*). From the expression of T in terms of C and since ||C|| = 1, it 
holds 

II^E||<II^II<II$H , < |i r -i||V2 

II dt II ~ II dt II ~ W^Wm II 

in virtue of the bound (|4.10p on K[$] using the fact that H'&Hl 2 = K. Inserting the 
last bound above in (|4.15|) and integrating over t yields 

lire*)- 1 1| < lir(o)- 1 1| + const. /Vw" 1 !! 372 ^, 

Jo 

for all t G [0,T*). Because of (|4~T2)) . this implies that f Q T * WT(s)- 1 1| 3 / 2 ds = +oo. 

So far we have proved the local well-posedness of the MCTDHF equations in 
X m for every m > 1 and the existence of a maximal solution in H 1 until time T* 
when the density matrix becomes singular. We prove now that T* is the maximal 
time of existence regardless of the imposed regularity on the solution. Let (C, <f>) 
be a solution in X2, then it is in particular a maximal solution in X\. We have to 
show that the H 2 norm of $ cannot explode at finite time < r < T* . Indeed, for 
any r < T* , we have 

(4.16) max ||r(t) -1 || < 1 

0<t<T 

by definition of T* , hence 

(4.17) max ||$(t)||m < 1. 

0<t<T 

From the Duhamcl formula for the PDEs system (|4.3p ~ (|4.4p and using the bounds 
flSH) and (TO]) together with ||$|| L 2 = 1 and ||C|| = 1, we get for all t G [0, r] 

||*(t)||fl* < \\<S>°\\ H i + C sup / ||$( S )|| ff , 

[0,r] JO 

where C is a positive constant that only depends on the local bounds (|4. 1 6(1 and 
(|4.17|) . By Gronwall's lemma we infer 

max \\<£>(t)\\ H 2 < e° T 

0<t<T 

hence the conclusion. The proof for any m > 2 follows then by a straightforward 
induction argument using the corresponding bounds (|4.6j) and ()4.8j) by assuming 
that max < t < T ||$(i)|| H m-i < 1. 

The proof of Theorem 14. II is now complete. 
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4.3. Existence of Standing wave solutions. In the present case the equations 
for the coefficients write (|3.19j) while Eqn. (|3.18|) for the orbitals becomes: 

(4.18) T(C) H $ + W[C, $] $ = A • $ 

according to Proposition 13.161 In [25] Le Bris has proved the existence of ground- 
states - that is, minima of the energy over the set Tn^k - for the physical Hamil- 
tonian on the whole space R 3 , and under the assumptions K = N + 2 and 

Ylm=i Zm > N — 1. Later on Friesecke extended this result to general admissible 
pairs (N,K), under the same assumption on the nuclear charge. Finally Lewin 
proved the existence of infinitely many critical points of the MCHF energy for any 
pairs (N, K) , hence the existence of infinitely many solutions to the coupled sys- 
tem (|4.18p - ()3.19|) that satisfy the full-rank assumption. All these solutions then 
give rise to infinitely many standing waves of the MCTDHF system and thereby to 
particular global- in-time solutions. 

The conservation of the invcrtibility of the matrix F(t) being an essential issue 
in the MCTDHF setting it is natural to give sufficient condition for such property. 

5. Sufficient condition for global-in-time existence 

In this section we focus again on the A^-body Schrodinger operator fll.lj) with 
physical interactions (|4.1I) . For any K > N + 1 with fixed N, we denote 

1(K) = inf : (C, $) G T N ,k) 

the "K- ground state energy" . Obviously one has 

(5.1) VK' < K < oo, inf a(H N ) < 1{K) < X(K'), 

with inf <t(Hn) being the bottom of the spectrum of Hn on L^(il N ). Recall that 
the maximal rank hypothesis corresponds to the following equivalent facts : 

(i) The rank of the operator [ir(C, <J>) (g> 7r(C, $)] : i is equal to K\ 

(ii) The K x K matrix F(C) is invertible; 

(iii) The smallest eigenvalue of T(C) is strictly positive. 

Since this is satisfied for K = N (Hartree-Fock case) and since K must be ad- 
missible, we now assume that K > N + 2. The main result of this section is the 
following: 

Theorem 5.1. Let (C7°,$°) G T NM be an initial data in ^T!\ with F(C°) in- 
vertible. Assume that T* < +oo then 

£(7r(C7°,$ )) >1{K - 1). 

As an immediate by-product we get a sufficient condition ensuring the global- 
in-time invertibility of the matrix IT(C(t)). 

Corollary 5.2. J/(C°,$°) G dT N . K satisfies 

(5.2) 1(K) < £ (7r(C7°, $ )) < I(K - 1), 
then T* = +oo; that is, the maximal solution is global-in-time. 

Remark 5.3. The hypothesis X)m=i z m > N in (|4.ip implies the relation I(K) < 
T{K — 2) |25U18| . Therefore can always be satisfied by changing K into K — 1. 

Remark 5.4. A key difficulty in the proof of above theorem is that the energy 
functional $ i— > £{^) is not weakly lower semi- continuous in H 1 (M. 3N ) while it is 
in H 1 ^ 31 *) for any bounded domain f2 as already observed by Friesecke |17j . When 
fl is a bounded domain of R 3 or when the potential U is non-negative, the proof of 
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Theorem \5.1\ is much easier thanks to the lower semi- continuity, and it is detailed 
in [3] ■ In the general case the proof is in the very spirit of Lewin's one for the 
convergence of critical points of the energy functional [26| . 

Proof of Theorem OjHJ Let (C, $) be the maximal solution to (pHT]) on [0, T*) 
with initial data (C°,<J>°) given by Theorem Ol We assume that T* < +00, then 



]imsup||r(C(t)) 

t/T* 



— 1 1 

= +OO. 



Equivalently, with the eigenvalues of r(C) being arranged in decreasing order < 
Ik < Ik-i < • ■ ■ < 7i < 1, this means 

liminf "fait) = 0. 

Then there exists a sequence t n converging to T* , a positive number (3 and an 
integer N + 1 < m < K such that 

(5.3) lim 7 m (£„)=0 and < /3 < 7 m _i(i„) . 

n — >+oo 

Indeed, since X)feLi7fe(*) = ^1 f° r au * e [0,T*) , at least N eigenvalues stay 
away from zero when t goes to T*. We denote C" = C(t„ ), $™ = $(£„), 7" = Ji(t n ), 
F n = T(C(t n )) and so on for other involved quantities. 

For all n > 1, (C™,$™) € 8J-n,k- Thus according to Proposition 12.51 there 
exists a unique sequence of unitary transforms U n G O^- that map (C n ,$ n ) into 
(C™, $'") with <£>'" being an eigenbasis for the operator 7" := 7*™. In particular 
the corresponding matrix F' n := T(C") is diagonal. In other words. 

*":=7r(C™ = = E c '^<" = 

(T a 
K K 

7" = E tf?®?J n = E ^ ® 

Since the group of unitary transforms is compact, we may argue equivalently on the 
sequence (C™, $' n ) that we keep denoting by (C",$ n ) for simplicity. From (|5.3j) 



(5.4) 



Then, 



lim 7" = for all m < i < K, 

n — >+oo 

liminf 7," > (3 > for all 1 < i < to- 1. 

n — '+00 



(5.5) for all cr e E#, {m, . . . , K] n a ^ =S> lim c£ = 0, 

Tl — * + 00 

for 7™ = X)ig,r I c ct| 2 m virtue of (|2.12p . In particular, the sequence C n £ S r ~ 1 
being compact 

(5.6) lim J2 K\ 2 = 1 - 

crc{l,...,m-l} 

We decompose 
with 



* n = 7r(C™, $ n ) = *+ + % 



*n = E C " $ - *« = E 

ern{m,...,A"}^0 o-n{m,...,A'}=0 

Then 

lim II*: " 



n— »+oo 



■ n llL 2 (R3™) 
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as a consequence of (|5.5[) and since each determinant <I>™ is normalized in L 2 (M. 3N ). 
Hence 



(5.7) 



lim llvp" - v& + H = f) 



Since the MCTDHF flow keeps the energy constant, we have 

£(7r(C n ,$ n )) =f(7r(C°,$ )), 

for all n > 1. This property provides with additional information on the sequence 
(C n ;<f> n ). Using the fact that the 0"'s diagonalizc 7™, the energy (|4.5|) rewrites 



if 

£(7r(C\ 7? / 



(5. 



if 
i=l 



|vci 2 + f/|CI 2 
|v^| 2 + c/|CI 2 



A" 



z,j,k,l=l 

dx, 



where in (|5.8|) we used the positivity of the two-body interaction potential v. By 
the Kato inequality, for any < e < 1, there exists C e > such that 



|C/| < -eA + C, 
in the sense of self-adjoint operators. Then 

K 

E7" 



Therefore, inserting into 



Et? / \VM\ 2 dx)-C £ N. 



K 

E^r / \V€\ 2 dx<cste. 



Thus, for all 1 < i < K , ^7™ >™ is bounded in ff^R 3 ). Then, from and 

extracting subsequences if necessary, we obtain the alternative 

(5.9) 

for all m < i < K, y/^yf ' <j>™ converges to weakly in iJ^R 3 ) and strongly in L 2 (M 3 ), 



and 

(5.10) for all 1 < i < m - 1, (f>? is bounded in H 1 ^ 3 ). 

Since, under the hypotheses on U, the map ip 1— > J R3 £/ |i^| 2 is weakly lower 

semi-continuous on iJ 1 (R 3 ), we deduce from (|5.9p that 

(5.11) 

- r l 

iv&t + y wr dx> 

n—*-\-oo 



> lim inf ]T 7 » / \\VW\ 2 +U\W 



lim inf 7™ 

n — >+oo ^ — ' 
»=1 

We now check that 
(5.12) 

A" ra-1 

hminf Yj 75H-D(^^;^^)=hmmf £ 1? jk i D \ <t>% 

i,j,k,l — l i.j,k,l—l 

by showing that 

(5.13) hminf £ 7 £ fei # ; ^ = 0. 



i , j , fc , Z = 1 

{i ) i,fe,i}n{m,... ) K}^0 
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Let {i,j,k,l} n {m,...,K} ^ 0. We assume without loss of completeness that 
i>m. From the expression (|2.10|) for 7 $.y, we observe that 



< 



mm I W 7 



7j J mm 



<min(7?;7^;7 J n ;7D 1/2 : 



(5.14) 

since < 7" < 1. We thus get 

(5.15) i£{i,j,k,l}n{m,...,K}^9, Urn 7 » fc ., 

n — '+00 J 

from (|5"1| . Then thanks to fl~T)) and ([5~ll|) 



< 



< 



0, 



7?V^I|V^|Ua||0?|U»||^|U»|I^IU' 

7? 



since the L 2 norms of the orbitals equal 1 and since in any case y/jj^cj)™ is bounded 
in L 2 independently of n. Therefore each term which appears in the sum in (|5.13p 
converges to as n goes to infinity thanks to (|5.9p . Claim (|5. 13[) then follows. 

Gathering together (|5.1ip and (|5.12|) we have 



liminf£(7r(C n ;$ n )) 

n — '+00 



m— 1 

(5.16)> liminf [ V 

i=l 



0, 



|vci 2 + ^ICP 



m— 1 
i,j,k,l— 1 

The point now consists in showing that the right-hand side in (|5.16[) is bounded 

Indeed, let us set $ n = 7r(C n ,$") where C n = 



from below by liminf Eify 

n — >+oo 



(cZ)«c{i,...,m-i} e CU-iJ and$" = (#,..., C-i) e O y{nr .. There is a slight 
difficulty arising here from the fact that (with obvious notation) 7™ is close but 
different from 7" i5^- and similarly for 7™- fci and 7 yjy. (Also C™ is not normalized 
in C m_1 (only asymptotically) but this will be dealt with afterwards.) 

First we observe that because of (|2.1ip for every i, j €{!,-■■ ,m— 1}, 



7" S n 



(<tUt) n{m,- ,K}^0 
iecrjer ,<?\{i}=r\{j} 



goes to as 11 goes to infinity thanks to (|5.5p . In addition, each term of the 



form / R3 |§V#- Vft + Uft-V} 
{1, • • • , m — 1}. Therefore 
(5.17) 



dx is bounded independently of n for z, j € 



m— 1 
i=l 



dx+o(l). 



For the same reason, and with obvious notation, for all 1 < k, I < m— 1, 

since according to (|2 . 1 0[) the extra terms in these differences only involve coefficients 
c™ with a n {m, . . . , X} ^ 0. Again each term of the form £>(</>" </>™ ; is 
bounded independently of n for fe, I 6 {1, . . . , m — 1}. Therefore 



(5.18) £ 75h^(^?^; 



m— 1 

i,j,k,l—l 



T3J w s(#0?;^#) + o(i). 
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Therefore, gathering together ([BTTB]) . |5~T7)) and ([5~T%|) . 

(5.19) liminf f (tt(C"; $")) > liminf S (ir(C n ; $")) = liminf £(*"). 

n — >+oo n — >+oo n — >+oo 

Since C" is not in S^m-i) -1 (it is only the case asymptotically), (C n ;$") is not 
in JFjv jTO _i, thus we cannot bound immediately £ (\& n ) from below by I(m — 1). 
Anyway, in virtue of (|5.6|) . 

(5.20) lim ||$ ,l || 2 = 1. 

n — >+oo 

The energy being quadratic with respect to "J 

(5.21) = ||^"|| 2 £(^^) > ||1<"|| 2 /(m- 1), 

^|| ^"||' 

for $'7||* n || G ^jv.m-i for all n > 1. Gathering together (^1^1 . (fFT^U|) and (|F^T|) 
and taking the limit as n goes to infinity we deduce 

(5.22) liminf SU{C n , $")) > I(m - 1). 

n — >+oo 

Hence the theorem. □ 

Remark 5.5. When fl is a bounded domain o/R 3 , any sequence in Fn,k is rela- 
tively compact in C r x L 2 (Q) K thanks to the Rellich theorem. On the other hand, 
the energy functional W i— > £ (\f r ) is weakly lower semi- continuous in if 1 ^ 3 ^). 
Therefore it is easily checked in that case that 

liminf £U{C n ;<S> n )) > liminf £(tt(C*; $*)) > I(m - 1) 

n — >+oo n — >-t-oo 

with (C*; $*) € TN,m-i being the weak limit of the sequence (C n ;$> n ) introduced 
in the above proof. 

Remark 5.6 (Stability, Consistency and Invertibility of the density matrix T). 
The main factor in the instability of the working equations or any gauge- equivalent 
system, is the inverse of the density matrix. In the present section, criteria for 
the global invertibility oflT(C) have been given. These criteria do not provide with 
an uniform estimate for ||T -1 ||, and furthermore increasing the consistency of the 
MCTDHF approximation leads to an increase of the number K of orbitals. As 
usual consistency and stability are both necessary and antinomic. Indeed, the most 
obvious observation is that one always has 

ii ii _ N , 

for IT has at most K positive eigenvalues whose sum equals N . Therefore the 
smallest can be at most N/K. These considerations lead either to a limitation on 
K or to a regularization or a "cut-off" of IT -1 . In fact the "consistency" in the 
sense of numerical analysis is obtained with fixed N by letting K go to infinity. 
This is basically different from the idea ( in spirit of statistical mechanics ) of letting 
N go to infinity [3] . 

6. Stabilization of T and existence of I? solutions 

In the above analysis, both for existence of maximal solutions and for global 
invertibility of the density matrix, the conservation of energy plays a crucial role. 
Besides the theoretical interest, the analysis of an MCTDHF system with infinite 
(or non conserved) energy but finite mass is relevant. Indeed, to circumvent the 
possible degeneracy of the density matrix, physicists resort to ad hoc methods like 
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perturbations of this matrix in order to ensure its invertibility. Typically, this is 
achieved as follows 

(6.1) T e =T + eId 
(see e.g. [7]), or by taking 

(6.2) r e = T + e cxp(-r/e) 

for small values of e (see [S]). Note that in latter case vanishing eigenvalues are 
perturbed at order e while the others arc unchanged up to exponentially small 
errors in terms of e. Then the perturbed system reads for an e > 

dC 



(6.3) 



1 5 = H * + iMC)- 1 (i - p*) w[c, 

at 

t C(0) = C°, $(0) = $°. 



On the other hand, when a laser field is turned on, the Hamiltonian of the system 
is then time-dependent which is a relevant configuration from the physical point of 
view (see [7J and Section [7J below). In such situation, the conservation of the energy 
fails and a recourse to alternative theories is necessary. 

However in both situations the L 2 norm (which corresponds to the electronic 
charge) is conserved and this justifies an L 2 analysis of the MCTDHF outside the 
energy space. Therefore the Strichartz estimates turn out to be a natural tool in 
the same spirit as Castella [§] and Zagatti In [21], existence and uniqueness 
of global-in-time mild solutions has been obtained for L 2 initial data. As in the 
previous section (and with the same notation) the perturbed working equations are 
written in "Duhamcl" form 



C{t)=C(0) + K[<£{s)}C{s)ds, 
Jo 

$(t) = S(t) $° - i [ S(t- s) U $(s) ds 
Jo 

-i [ S(t- s)T t [C(s)}~ 1 (I-P*)W[C,*]$ds, 

where S(t) = exp[—^itA] denotes the group of isometries generated by — | A on 
L 2 (K 3 ,C). The potentials U and v = v{\x\) belong to hi + L°° . 
From the relation 

and 

\\S{t)<l>\\L*(w?) = HWL^m) 
one deduce by interpolation the so-called Strichartz estimates 

||5(t)0||iP(o,T;i9(K3)) < C(q)Ti 2 , 

that hold for any Strichartz pairs (p,q) G [2, +oo] x [2,6] with ^ = (5 — |)- 
(Strichartz estimates for the endpoints p = 2 and q = 6 are more intricate and due 
to Keel and Tao [22]). 

Following Zagatti (35] and Castella [H], the spaces 

X T = L°°(0,T;C r ) x (L oo (0,T;L 2 (R 3 ))nU'(0,T;L g {R 3 )) Xl 
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are introduced for any Strichartz pairs. For some R > and some T > small 
enough, the non-linear operator (C, <!>) i— ► L(C, $) which appears in the Duhamcl 
integral 



mm) 



it 1 \ 

*'s)] C(s)ds 



V 



J S{t - s) (u $(s) + r e (a) _1 (J - P#) W[C(s), $(s)] $(s)) ds 



is a strict contraction in the ball 

{(C, $) G X T : ||C||c + ||$||l°°(0,T;L 2 (R 3 )) + \\®\\lp(0,T;Li(R 3 )) < -R} ■ 

Next using the conservation of the L? norms of the orbitals and the estimate 

II^IIl^o.til^r 3 )) iS II^°IIl 2 (r 3 ) 

one follows the lines of Tsutsumi in [33] to get existence and uniqueness of a strong 
solution in (see the details in |31j). This is summarized in the 

Proposition 6.1. Let e > 0. For any initial data (Co,$o) G OJ-n^k and for any 
Strichartz pairs (p,q), the e-regularized working equations admit a unique strong 
solution 

(C e (t),* e (t)) G L°°(R + ;C r ) x (L°°{R + ;L 2 {M 3 ))nLf oc (R + ;L q {R 3 )) K 

that lives in Tn.k for all t > 0. If in addition <& G H 1 (R 3 ) K then $ e (t) G 
C ,0 (R+; J ff 1 ( R3 )) if - 

Eventually one expects that whenever the original solution is well-defined (with 
a non degenerate density matrix IT(t)) on a time interval < t < T* it will be 
on the same interval the limit for e — *■ of the solution of the perturbed working 
equations. This is the object of the following 

Theorem 6.2. Let (C ,$o) G S^ 1 x (H 1 ^ 3 )) 1 *. -4ssMme the corresponding 
solution (C(t),<&(t)) to (|3.27|) is well-defined on [0, T] and zs suc/i £/ia£ 

(6.4) sup ||r(t) _1 | < A/ < +oo. 

0<t<T 

Then, on the same time interval it is the limit in C r x H 1 (M. 3 ) K for e — > o/ £/ie 
solution (C e ,$ e ) £o t/ie regularized problem (|6.3[) mi/i same initial data. 

Proof. We first recall the obvious a posteriori bounds 

yen = yen = i, ii$!| i2 = n$ e || L2 = i 

on [0, T], and, as a consequence of (|6.4p and the energy conservation, 

max \\$(t)\\ H i < M' 

0<t<T 

with M' = M'(£(ir(C ,<f> )),M). We can also rely on the orthonormality of the 
orbitals in <£> and <& e . We introduce the notation 

*.-(%). ^U). «w>-(®> 



where 

= (I-P*)W[C,*]$ 
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and where the index (e) means that the claim holds both for the regularized system 
and the initial one, uniformly in e. System (|6.3|) can also be written in synthetic 
form: 

(6.5) U(t) = e- itA U - i f e- l ^-^ A B{U{s)) ds, 



(6.6) U e (t) = e~ itA Uo - i / e-^- s ^ A B e (U e {s)) ds 



Since the initial < & e (0) = <I>o is in H 1 and since the regularized system propagates 
the regularity, <!> e is in H 1 (M 3 ) K for all time. 

We fix e > 0. We introduce a parameter 77 > to be made precise later and the 

set 

7 e = {iG[0,T] : \\U e (t)-U(t)\\ <!!} 
with \\U\\ = 1 1 C 1 1 + The mappings t >—> U/ e \(t) being continuous from [0,T] 

toX T := C r x L°° (0,T; H 1 ^ 3 )*) , the set I e is closed and since ||J7 e (0)-L/(0)|| = 0, 
there exists a maximal time T t > in I t such that 

Vie [0,1;], \\U e (t)-U(t)\\ <r). 
We now prove by contradiction that T c = T. Assume then T e < T. 
Subtracting (|6.5|) to (|6.6p and taking norms first yields to 

\\C e (t) - C(t)\\ < f ||C e ||||K[$ £ ]-K[$]|| + ||C e -C||||K[$]j|, 
Jo 

(6.7) <C( V ) [ \\U £ -U\\ds 

Jo 

for all < t < T e . Here and below C(-q) = C(M,£ (tt(C , $0)), v) denotes a positive 
constant that may vary from line to line but that is independent of e and continuous 
and non-decreasing with respect to 77. Indeed we use the fact that the non-linearity 
$ 1 > K[$] is locally Lipschitz continuous in H 1 (Subsection 14. ip together with the 
uniform bound 

max ||$ellffi < M' + n. 

0<t<T e 

On the other hand, we write 

ll*e(*) - *(*)ll < f WnuMm \\T~ e l T- l \\ + lir^H \\M(U £ ) M(U)\\m 
Jo 

(6.8) < C{rj) j (HIT7 1 - T- 1 !! +\\V e - U\\ H1 ) ds, 

by using the local Lipschitz bounds of U 1— > M(U) given in Subsection 14.11 We 
now turn to the quantity 1 1 IT^ 1 — T -1 ||. Both regularization (|6.1|) and ()6.2|) of the 
density matrix take the form : 

r e -r(c e ) + e 5 (c £ ) 

with ||fif(C e )|| < e. Then, 

||r e -r|| < ||r(c e )-r(c)|| + e 

(6.9) <K(\\C e -C\\+e) 

by using the obvious bound ||r(C)|| < ||C|| 2 for C, C e G S r ~ 1 , where k only depends 
on iV and K. We now assume that 

1 



(6.10) e, 77 < 



4kM' 



MCTDHF EQUATIONS 40 

where M is given in the statement of the theorem. Using 

r e = (i- (r - r £ ) r^ 1 ) r, 

we deduce 

17 1 = r~ 1 (1- (r-r^r- 1 ) 1 = r~ 1 ^ ((r-r^r- 1 )" 

n>0 

Therefore 

H r -i _ r -i|| < ^ \\f - r € \\ n yr- 1 !!^ 1 

n>l 

< ^M n+ V l (||C e -C|| +e) n 

n>l 

by using (|6.9p . Hence 

|| r -i _ r -i|| < M 2 K _ C || + e ) ^ M n K n Q| Ce _ C || + £) „ 

n>0 

(6.11) < 2M 2 n(\\C e -C\\ +e) 

since M k (\\C e (t) - C(t)\\ + e) < \ by and for t in [0, T e ). Inserting ([(HI]) in 

we get: 



(6.12) \\^{t)-m\\<C(v) / (\\U e (s)-U(s)\\+e)ds. 

Jo 

Eqn. ([6~T2]) together with (pT7)) finally leads to 

(6.13) \\U e (t) - U(t)\\ < C{r,) I (\\U e (s) - U(s)\\ + e) da, 

Jo 

for all t £ [0, T c ]. Eventually thanks to Gronwall's inequality, 

(6.14) max \\UJt) - U(t)\\ < e e CM T . 

0<t<T e 



With r) as in (|6.10[) . next 

p <T min I — 

MkM' 2 



(6.15) e < min( 1 ^ e -c(,)T ); 



we get 

max \\UJt) - U(t)\\ < T j-. 

By continuity of 1 1— ► || J7 e (t) — J7(t) || , we may then find T e ' > T e such that max || Z7 e (i) - 

C7"(i) || < r). Hence the contradiction with the definition of T e . Therefore, I e = [0, T] 
and, going back to (|6.14[) we obtain: 

(6.16) max \\UJt) - U(t)\\ < ee c ^ T , 

0<t<T 

for say r\ = 4k 1 a/ and e small enough, satisfying (|6.15[) . whence the result. □ 

In the forthcoming (and last) section we comment on straight extensions of the 
above analysis. 
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7. Extensions 

The present contribution is focused on the algebraic and functional analysis prop- 
erties of the MCTDHF equations for fcrmions. Multi-configuration approximations 
can also be considered for symmetric wave- functions or also for wave- functions with 
no symmetry (see e.g. [5j [23] ) - The mathematical analysis of the equations which 
play the role of the "working equations" of Section [3] is similar. On the other hand, 
the fermionic case is important by itself and leads to much better geometric struc- 
ture in terms of principal fiber bundle as described in Section [2] Hence our choice. 
Our results could be generalized to general (symmetric) n-body interactions as well 
including the n-body density matrices. 

7.1. Beyond Coulomb potentials. Although above results and proofs are mainly 
detailed for Coulomb potentials they carry through more general real- valued poten- 
tials. Indeed well-posedness results in H 1 and H 2 are still valid for U and v in the 
class L p (R 3 ) +L QO (R 3 ) with p > 3/2, and v > 0. These conditions ensure that H N 
is self-adjoint in L 2 (Q N ), that the one-body operator — |A + U is a semi-bounded 
self-adjoint in L' 2 (M. 3 ) with domain H 2 (M. 3 ) and that the Kato inequality holds for 
the potential U. Under these assumptions, the energy space is C x H 1 (M 3 ) K (re- 
spectively C r x H£{n) K when fi is a bounded domain) and the propagator e 

is a one-parameter group of unitary operators in H 2 (M. 3 ) and in H 1 (M 3 ). 

For the global well-posedness sufficient condition to hold true (Theorem 15.11 and 
its corollary) further conditions on the potentials are required to ensure that the 
energy functional is weakly lower semi-continuous on the energy space. Sufficient 
conditions are (for example) U > or U— (the negative part of U) tending to at 
infinity at least in a weak sense. 

7.2. Extension to time-dependent potentials. One of the basic use of the 
MCTDHF is the simulation of ultra-short light pulses with matter [37] . Describing 
this situation leads to the same type of equations but with the one-body Hamilton- 
ian H being replaced by a one-body time dependent hamiltonian 

H^aW := (*V + A(t)) 2 + w(t) U(x) 

with oj(t) and A(t) real, A(Q) = and U as in the above subsection. A typical 
example is A(t) = Ao exp ( — (t/r) 2 ) sin(crf) for some positive real parameters Aq, 
a and r (37J [38]. This does not change neither the algebraic and geometrical 
structure of the equations nor the definition of the density matrix IT nor the notion 
of full-rank. The potential vector A being independent of the x variable the energy 
space is H . With convenient hypotheses (say u> and A continuous, bounded with 
bounded derivatives), the results in Section H] concerning local- in-time H 1 well- 
posedness of the Cauchy problem remain valid. For generalization of the use of 
Strichartz estimates and the local L 2 well-posedness one should follow for example 
|10| . Since the energy is now time-dependent extra hypothesis have to be introduced 
for the persistence of the full-rank assumption done in Section [5] 

Assume that ui(t) and A(t) take their values in a bounded set (the set of "control" 
C) and that their derivatives are also bounded. The system So Q3.25P with H 
replaced by H u keeps on preserving the constraints since Lemma [3~2l only relies 
on the self-adjointncss of the Hamiltonian. Similarly solutions to (|3.25p satisfy the 
Dirac-Frenkel variational principle. The energy is no longer conserved by the flow. 
Indeed, following the lines of the proof of Corollarv l3.41 we have 

±£(*(t)) = ±(H(t)*(t)\m) = (u/(t) + 2A(t)A'(t)) <*(t)|*(t)>, 
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with the prime denoting time derivatives. However 

£(*(t)) =u(t) + 2A(t) 2 +£(V(0)) -u(0), 

and the energy in controlled for any finite time, whence the existence of a maximal 
solutions in H 1 as long as the matrix F(C(i)) remains invertible. 

To adapt the result concerning the global full-rank hypothesis, we introduce the 
minimization problems for any real numbers Q and A 

with 

^,i(*) = ((H^r + lw[c, 

\ V Z 7 J L 2 (Q.) K 

for ^ = 7r(C, $). The global-in-time conservation of full-rank in Theorem 15.11 
remains true under the hypothesis 

£ U (0),A(0) (*(0)) <m£{X 0iA (K-l) : |o)| < |MU-(r +) , |A| < ||A|| £ co (B+) } 
- ||w|| £ cc( R+) -2||A|| ioo(R+) +w(0). 

There is a lot of room for improvement in the above argument. For example, if we 
assume that, for all time, the solution \& = 7r(C, $) € OBn.k satisfies 

<^*(*)|*(t)> </*(*) <tt(t)*(*)|*(t)) 
for a given function /i, then by the Gronwall lemma 

(«(t)#(t)|tf(t)) - (W(0)*o|*o> ^ cx p(/ 
The result of Theorem 15.11 remains true provided 

£(tfo) = (W(0)* |*o> < I(A' - 1) - cxp ( £ h(s)d S y 

7.3. Discrete systems. The emphasis has been but in particular for the functional 
analysis on the case when f2 = R 3 although in the first part we have described the 
problem in any open subset of M 3 . In fact all the formal and algebraic derivations 
can also be adapted to the case when Q is a discrete set equipped with a discrete 
Lcbesgue measure and in particular when Q is a finite set Such situation is 
important for two reasons. On the one hand many models of quantum physics 
(the Ising model for instance) involve a discrete Hamiltonian defined on a discrete 
set. On the other hand the discretization of the original problem in view of any 
numerical algorithm leads to a discrete problem. 

Up to now only a rough a posteriori error estimate has been proven. However 
if the MCTDHF algorithm is applied to a discrete model say of dimension L then 
one always has K < L. The error formula (|3.13|) shows that, for K = L, the 
MCTDHF algorithm is exact. It should be eventually observed that in general 
the two operations : - Discretization of the original iV-particlc problem and use 
of a MCTDHF approximation or - Use of a MCTDHF approximation and then 
discretization of the equations, lead to different algorithms. 

Appendix - Proofs of technical lemmas in Subsection 13.31 

Proofs of Corollary 13.71 and Lemma 13.111 For a and r given and fixed 1 < 
j < N it is convenient to denote by l] a T ^ the column vector in C N with entries 

(^<r(i),r{j))l<i<N andb y 

[U T (i),U T (2), . . . ,U t (jv)]<t 
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the determinant composed with these vectors. With this notation (13. 14)) gives 

(7-1) ^^ = X>, TW U CT , fe 

fe=i 

Differentiating the relation 

Ucr, r = [U T (X),U T (2), ■ ■ ■ , C^r(JV)]cr 

and using the multi-linearity with respect to the column vectors and Eqn. (|7.ip one 
obtains: 

(7.2) i ^ T = Mk tT (j)[U T (x), U T (2), . . . , U T (j-i), Uk, J7 T (j+i), . . . , f7 r (jv)]o-- 

i<fc</f 

l<j<N 

On the other hand since U(t) is a flow of unitary matrices it is solution to a differ- 
ential equation of the following type: 

(7.3) i^- =J2i U T'(i),U T , (2) ,...,U T , {N) ] a M T , <T 

r' 

Identification of the coefficients of 

\U T i(x), U T i(2), ■ ■ ■ , U T i( N )] a 

gives, taking in account the number of permutation needed to change 

r(l), r(2), . . . t(j - l),k, t(j + 1), . . . r(N) into r'(l), r'(2) . . . , r'(N) 



(A') 



ker'jer 
r'\{fe}=r\{j} 



Let us now prove (|3.2ip . Let a, r G *£>n,k- We first observe that 

JV N 



(7.4) ^G^^^^A-.-AGij,, 



A ... A 



>a{N)- 



Now we use (|2.5[) and the Laplace method to develop a determinant with respect 
to the row that contains the terms involving G to get 

JV N 

5^(G X4 $ ff |$ T ) =$^(&r(i) A...AG^ (j) A...A^ (JV) |$ T ) 

i=l i=l 
A 

= J! ( G ^<7(i).^r(j)> ^\{a(i)},r\{T(j)}) 

in virtue of (|2l| . Hence (j3~2l"j) using (|5TTS|) and the definition of M. □ 

Proofs of Theorem l3~13l and Theorem Let (C(i),$(t)) be a solution 

to iSo and let G be as in the statement of the theorem. With My = (G (f>i,(j}j) 
we define the family of unitary transforms U(t) according to Lemma 13.111 and 
d(U)(t) = U(i) is then given by Corollary O We set V = U, C'(t) = V(t) C(t) 
and $'(t) = U(t)<Z>(t). Thanks to ([335]) . V solves 

r rfv — 

i — = -VM. 
(7.5) { dt 



[v(0) = d(Z7°). 
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Then, for all a £ Hn,k, 

i ^- = i^-C + Vi^- = -VMV* C' + V * I V c *) 
at at at x 1 ' 

= -VMV* C +Y(H^ I V C '*V) 
thanks to (|2.25p and (|7.5[) . On the one hand, since V is unitary, 

On the other hand, when M is obtained through G, we get by a direct calculation 
from (pHI]) 



N N 
r i—l i— 1 



= E(E G - < Wk = (E G - VI ' 



Combining these two facts we get the first equation in iSq, namely 
i — = (n *| Vc*) - ( E G ^ *| Vc* 



We turn now to the equation satisfied by To simplify the notation we use the 
shorthand T for T(C) and T' for T(C") respectively. Then, using T' = U F V* 
and (|3.14|) . we have 

a$' rfr/ a<3> 

at at at 

9$ 

= T' UMU* & + UFi — 
at 

= r' umu* $' + (i — p*0 [/ v*** [w *] 

(7.6) = r' UMU* $' + (I - P*,) V*/¥* [w *] 

thanks to ([2~26|) and since clearly P$< = P$ for Span{<f>} = Span{<3?'}. It is easily 
checked that when M is given through G we have 

(UMU*) ij = (G^,4>' j ) 

and therefore 

UMU* = P$/G$'. 

Hence (17.61) also writes 



9$' 

ir' — = r' g$' + (i-p*/)v*/**[w*] -(i-p*/)r'G$'. 



at 

We now check that, for all 1 < i < N, 

N 



a** 

A 7 - 

1 j=i 

thereby proving that 

N 

- = F' CI 4/ -4- (f - T>,.A Pw v& - 

at 



(i-p„)(rG$') i = (i-p $ ,)^ [E G «** 



a$' r 
i T' — = JT' G + (I - P«/) V$'** W * - J] G^ * 
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Indeed, for all £ £ L 2 (Q), using (|2~2"g)) in Lemma EE in ([77)1 and using fTi]) in 
(|7.8p . we have 

if 

<(i-p«) (r g $).,£) = ^ r ifc (G <mi -p*K) 



fc=i 

A" 



(7.7) =E(tuT[ G ^I^-[( I - P *)^ 



it=i - Pk 



N 



(7-8) = (E G -i 

-<d-P.)^[gG. J .], { 

by the definition (|2.24[) of §§: ', whence the result since £ is arbitrary in L 2 (S1). □ 
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